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The two-time Green function method in quantum electrodynamics of high-Z few-electron atoms 
is described in detail. This method provides a simple procedure for deriving formulas for the energy 
shift of a single level and for the energies and wave functions of degenerate and quasi-degenerate 
states. It also allows one to derive formulas for the transition and scattering amplitudes. Application 
of the method to resonance scattering processes yields a systematic theory for the spectral line 
shape. The practical ability of the method is demonstrated by deriving formulas for the QED 
and interelectronic-interaction corrections to energy levels and transition and scattering amplitudes 
in one-, two-, and three-electron atoms. Numerical calculations of the Lamb shift, the hyperfine 
splitting, the bound-electron g factor, and the radiative recombination cross section in heavy ions 
are also reviewed. 
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I. INTRODUCTION 

A great progress in experimental investigations of high-Z few-electron systems (see, e.g., [0,0) stimulated theorists to 
perform accurate calculations for these systems in the framework of quantum electrodynamics (QED). The calculations 
of QED and interelectronic-interaction corrections in high-Z few-electron systems are conveniently divided into two 
stages. The first stage consists in deriving formal expressions for these corrections from the first principles of QED. 
The second comprises numerical evaluations of these expressions. The present paper is mainly focused on the first 
stage. As to the numerical calculations, we give only a short overview of them in this paper. For more details we 
refer to §-f|. 

Historically, the first method suitable for deriving the formal expressions for the energy shift of a bound-state 
level was formulated by Gell-Mann, Low, and Sucher ||J|. This method is based on introducing an adiabatically 
damped factor, cxp (— A|£|), in the interaction Hamiltonian and expressing the energy shift in terms of so-called 
adiabatic S\ matrix elements. Due to its simple formulation, the Gell-Mann-Low-Sucher formula for the energy 
shift gained wide spreading in the literature related to high-Z few-electron systems pQ-O]. However, the practical 
use of this method showed that it has several serious drawbacks. One of them is the very complicated derivation of 
the formal expressions for so-called reducible diagrams. By "reducible diagrams" we denote those diagrams where an 
intermediate-state energy of the atom coincides with the reference-state energy. (This terminology is quite natural 
since it can be considered as an extension of the definitions introduced by Dyson |0| and by Bethc and Salpeter 
Jl9| to high-Z few-electron atoms.) As to irreducuble diagrams, i.e. those diagrams where the intermediate-state 
energies differ from the reference-state energy, the derivation of the formal expressions can easily be reduced to 
the "usual (A = 0) S'-matrix" elements in each method, including the Gell-Mann-Low-Sucher method as well (see, 
e.g., |12|,|l6|). Another serious drawback of the Gell-Mann-Low-Sucher method is the fact that this method requires 
special investigation of the renormalization procedure since the adibatic S^-matrix suffers from ultraviolet divergences. 
The adiabatically damped factor, exp(— A|i|), is non-covariant and, therefore, the ultraviolet divergences can not be 
removed from S\ if A ^ 0. However, from the physical point of view one may expect these divergenes to cancel each 
other in the expression for the energy shift. Therefore, they may be disregarded in the calculation of the energy 
shift for a single level. For the case of degenerate levels, however, this problem remains since we can not expect 
that the standard renormalization procedure makes the secular operator finite in the ultraviolet limit |0,0]. In 
addition, we should note that at present there is no formalism based on the Gell-Mann-Low-Sucher approach which 
would provide a proper treatment of quasidegenerate levels. Also no formalism in the framework of this approach 
was developed for calculation of the transition or scattering amplitudes up to now. The same difficulties emerge in 
the evolution operator method developed in Refs. EQ-pJ]. Attempts to solve some of these problems by modifying 
the Gell-Mann-Low-Sucher method were recently undertaken in |25|-|27|| . 

Another way to formulate a perturbation theory for high-Z few-electron systems consists in using Green functions. 
These functions contain the complete information about the energy levels and the transition and scattering amplitudes. 
In this way, the renormalization problem does not appear, since Green functions can be renormalized from the very 
beginning (see, e.g., ^8[). Up to now, various versions of the Green function formalism were developed which differ 
from each other by the methods of extracting the physical information from Green functions, i.e. the energy levels 
and the transition and scattering amplitudes. One of these methods was worked out in 29 -33). It was successfully 



employed in many practical calculations p3E6f. Since one of the key elements of this method consists in using 
two-time Green functions, in what follows we will call it the two-time Green function (TTGF) method. This method, 
which provides a solution of all the problems appearing in the other formalisms indicated above, will be considered 
in detail in the present work. 

As to other versions of the Green function method 1 13 17| , |47| -p5| , a detailed discussion of them would be beyond the 
scope of the present paper. We note only that some of these methods are also based on employing two-time Green 
functions but yield other forms of perturbation theory. In Refs. |13] , f4q -[5l| , the two-time Green functions were used for 
constructing quasipotential equations for high-Z few-electron systems. This corresponds to the perturbation theory in 
the Brillouin-Wigner form while the method of Refs. fl29|-[33fl yields the perturbation theory in the Rayleigh-Schrodingcr 
form. Various versions of the Bethe-Salpeter equation derived from the 2A^-time Green function formalism for high-Z 
few-electron systems can be found in P3|,p2| . In J5^j54|] the perturbation theory in the Rayleigh-Schrodinger form is 
constructed for the case of a one-electron system where the problem of relative time coordinates for the electrons does 
not occur. 

The relativistic unit system ( % = c = 1 ) and the Heaviside charge unit (a = f-,e < 0) are used in the paper. 



II. ENERGY LEVELS OF ATOMIC SYSTEMS 

In this section we formulate the perturbation theory for the calculation of the energy levels in high-Z few-electron 
atoms. In these systems the number of electrons denoted by N is much smaller than the nuclear charge number 
Z. It follows that the interaction of the electrons with each other and with the quantized electromagnetic field is 
much smaller (by factors 1/Z and a, respectively) than the interaction of the electrons with the Coulomb field of the 
nucleus. Therefore, it is natural to assume that in zeroth approximation the electrons interact only with the Coulomb 
field of the nucleus and obey the Dirac equation 

(-ia • V + 0m + Vc(x))V„(x) = e n ip n (x) . (2.1) 

The interaction of the electrons with each other and with the quantized electromagnetic field is accounted for by 
perturbation theory. In this way we obtain quantum electrodynamics in the Furry picture. It should be noted that 
we could start also with the Dirac equation with an effective potential T4s(x) which approximately describes the 
interaction with the other electrons. In this case the interaction with the potential SVfe) — Vc(x) — V c s(x) must be 
accounted for perturbatively. Using the effective potential provides an extension of the theory to many-electron atoms 
where, for instance, a local version of the Hartree-Fock potential can be used as I4ff(x). However, for simplicity, in 
what follows we will assume that in zeroth approximation the electrons interact only with the Coulomb field of the 
nucleus. 

In the present paper we will mainly consider the perturbation theory with the standard QED vacuum. The transition 
to the formalism in which closed shells are regarded as belonging to the vacuum is realized by replacing iO with — iO 
in the electron propagator denominators corresponding to the closed shells. 

Before formulating the perturbation theory for calculations of the interelectronic interaction and radiative correc- 
tions to the energy levels, we consider standard equations of the Green function approach in quantum electrodynamics. 

A. 2iV-time Green function 

It can be shown that the complete information about the energy levels of an iV-electron atom is contained in the 
Green function defined as 

G(x[, ...x' N ; Xl ,... x N ) = (0|2V(&i) • • • ^(x' N )^(x N ) ■ • >(xi)|0) , (2.2) 

where ip( x ) is the electron-positron field operator in the Heisenberg representation, ip( x ) — i> T°> an d T is the 
time-ordered product operator. The basic equ ations of quantum electrodynamics in the Heisenberg representation 



are summarized in Appendix A. Equation (2.2) presents a standard definition of the 2A-time Green function which 
is a fundamental object of quantum electrodynamics. It can be shown (see, e.g., [|28||56|]) that in the interaction 
representation the Green function is given by 



G(x' 1 ,. ..x' N ;xi,.. .x N ) 

(Q\Tip- m (xi) ■ ■ -Tpinix'^ip^XN)- ■ ■ ^ in (xi) exp {-i J d 4 z Hi(z)}\0) 



(2.3) 



(0|Texp{-jfd 4 zft/Cz)}|0) 

{ Z^2 m \ I d4yi ' ' ' d4ym (°\ T ^( X 'l) • • • ^in{x' N )^ in (x N ) ■ ■ ■ 1p in (xi) 
m=0 ' •* 

'Mi{yi) ■ ■ ■Hi(y m )\0)}{j2 ^p / rf4z i ' ••<***! (0\THi(z x ) ■ ■■H I (z l )\Q)} (2.4) 
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ni(x) = - ® ta (xft M , An{x)]At n {x) - — [iMx), <M*)1 (2-5) 



where 



is the interaction Hamiltonian. The commutators in equation ( |2.5[ ) refer to operators only. The first term in ( |2.5| ) 
describes the interaction of the electron-positron field with the quantized electromagnetic field and the second one 
is the mass renormalization counterterm. We consider here that the interaction of the electrons with the Coulomb 
field of the nucleus is included in the unperturbed Hamiltonian, i.e. the Furry picture. However, there is also an 
alternative method to get the Furry picture. In that method the interaction with the Coulomb field of the nucleus is 



included in the interaction Hamiltonion and the Furry picture is obtained by summing infinite sequences of Feynman 
diagrams describing the interaction of the electrons with the Coulomb potential. As a result of this summation, the 
free-electron propagators are replaced by bound-electron propagators. This method is very convenient for studying 
processes involving continuum-electron states. It will be used in the section concerning the radiative recombination 
process. 



The Green function G is constructed by perturbation theory according to equation (2.4). This is carried out with 
the aid of the Wick theorem (see, e.g., J28|). According to this theorem the time-ordered product of field operators is 
equal to the sum of normal-ordered products with all possible contractions between the operators 

T(ABCD ■■■) = N(ABCD ■ ■ ■) + N(A a B a CD ■ ■ ■) + N(A a BC a D ■ ■ ■) 

+ all possible contractions, (2-6) 

where N is the normal-ordered product operator and the superscripts denote the contraction between the correspond- 
ing operators. The contraction between neighbouring operators is defined by 

A a B a = T ( AB ^ _ tv(AB) = (0|T(AB)|0) . (2.7) 

If the contracted operators are boson operators, they can be put one next to another. If the contracted operators are 
fermion operators, they also can be put one next to another but in this case the expression must be multiplied with 
the parity of the permutation of the fermion operators. Since in the Green function the vacuum expectation value is 



calculated, only the term with all operators contracted remains on the right-hand side of equation (2.6). In contrast 
to the free-electron QED, in the Furry picture the time-ordered product of two fermion operators must be defined 
also for the equal-time case to obtain the correct vacuum polarization terms. As was noticed in |57|, the definition 

T[A(t)B{t)] = l -A{t)B{t) l -B{i)A{t) (2.8) 

provides the following simple rule for dealing with the interaction operator. It can be written as 

Hi(x) = eip- m (x)7 IJ ,ipi n (x)A% n (x) - Smtfj in (x)i}i n (x) (2.9) 

and then the Wick theorem is applied with contractions between all operators, including equal-time operators. We 
note that the problem of the definition of the time-ordered product of fermion operators at equal times does not 
appear at all if the alternative method for obtaining the Furry picture discussed above is employed. 

The contractions between the electron-positron fields and between the photon fields lead to the following propagators 

(0|T^ n (a# in (y)|0) = ± [°° du E ^f^ expHa;^-;/ )] (2.10) 



and 



^XMi^m= ■<•/' /70 cxp[ ^T y)] - (2-ii) 



Here the Feynman gauge is considered. In equation ( 2.10]) the index n runs over all bound and continuum states. 



The denominator in equation ( |2.3| ) describes unobservable vacuum-vacuum transitions and, as can be shown (see, 
e.g., p8[), it cancels disconnected vacuum- vacuum subdiagrams from the numerator. Therefore, we can simply omit 
all diagrams containing disconnected vacuum- vacuum subdiagrams in the numerator and replace the denominator by 
1. 

In practical calculations of the Green function it is convenient to work with the Fourier transform with respect to 
time variables, 

G((p?\ x' x ), . . . , (p'l x' N ); (p°, Xl ), . . . , (p° N , XJV )) 

/CO 
dxl ■ ■ ■ dx° N dxf ■■■dx% 
-OC 

x exp (ipfxf + ■■■ + ip%x% - ip\x\ ip%x° N ) 

xG(x' x ,...,x' N ;xx,...,x N ) . (2.12) 

For the Green function G((p' 1 °, x.[), . . . , (p%, x-' N ); (pi, xi), . . . , (p%, Xj\r)), the following Feynman rules can be derived: 
(1) External electron line 



y £?S(w, x, y) 



where 



S{u,x,y) =J2 



-0n(x) are solutions of the Dirac equation (2.1) 
(2) Internal electron line 



W -£n(l - *0) 



(2.13) 



(3) Disconnected electron line 



(4) Internal photon line 



x y ^j-oo^^v^-y)' 

where, for zero photon mass, D pa (uj, x — y) is given by 



D pCT (w,x-y) 



-9 P o 



dk exp (ik • (x — y)) 



(2tt) 3 w 2 - k 2 + (0 



in the Feynman gauge and by 

Doo(w,x-y) = 
A/(w,x-y) = 



1 



Ao = Ak = (1 = 1,2,3), 

(< 



47r|x — y 

dk exp(zk- (x-y)) / hh . 

(2^ - 2 -k 2 +,o r<-i^ m= l > 2 ' 3 ) 



(2.14) 

(2.15) 
(2.16) 



in the Coulomb gauge. (In this work we assume that the Coulomb gauge is used only for diagrams which do not 
involve a renormalization procedure. The renormalization in the Coulomb gauge is discussed in Refs. |58|j59|] .) 
(5) Vertex 




— 2itie^ p 8(bJi — LU2 — W3) / dx . 
(6) The mass counterterm 



X /\ 



2iri5(u> — uj')Sm J dx . 



(7) Symmetry factor (— l) p , where P is the parity of the permutation of the final electron coordinates with respect 
to the initial ones. 

(8) Factor (—1) for every closed electron loop. 

(9) If, in addition, an external potential <5V(x) is present, an additional vertex appears, 



-2irij 5(uj - u') J dx5V(x) . 

In principle, the Green function G contains the complete information about the energy levels of the atomic system. 
This can be shown by deriving the spectral representation for G. However, it is a hard task to extract this information 
directly from G since it depends on 2{N — 1) relative times (energies) in the time (energy) representation. As we will 
see in the next section, the two-time Green function defined as 

G(t', t) = G(t[ = t' 2 = ■ ■ ■ t' N = t'; t x = t 2 = ■ ■ ■ t N = t) (2.17) 

also contains the complete information about the energy levels, and it is a much simpler task to extract the energy 
levels from G. 

B. Two-time Green function (TTGF) and its analytical properties 

Let us introduce the Fourier transform of the two-time Green function by 

1 1 r°° 

Q(E]x\,...x! N \x. 1 ,...x N )5(E-E') = — — / dx°dx'° exp (iE'x'° - iEx°) 

2m N\ J_ 0o 

x<0|2>(x'°,xi) • • ■^(x m , X ' N Wx°,x N ) ■ ■•^(z°,x 1 )|0) , (2.18) 

where , as in fl2.2p, the Heisenberg representation for the electron-positron field operators is used. Defined by equation 



(2.18) for real E, the Green function Q can be continued analytically to the complex E plane. Analytical properties 
of this type of Green functions in the complex E plane were studied in various fields of physics (see, e.g., [|60|-|62|). 
In quantum field theory they were considered in detail by Logunov and Tavkhelidze in p3[ (see also J64[), where the 
two-time Green function was employed for constructing a quasipotential equation. To study the analytical properties 
of the two-time Green function we derive the spectral representation for Q . Using the time-shift transformation rule 
for the Heisenberg operators (see Appendix A) 

^{x°, x) = exp {iHx a )i;(0, x) exp (-iHx a ) (2.19) 

and the equations 

H\n)=E n \n), ^|n)H=^, (2-20) 

n 

where H is the Hamiltonian of the system in the Heisenberg representation, we find 

g{E;x[,...,x' N ;K U ...,x N )S(E - E') 

1 1 I" 00 
= 7T- T7T / dx°dx ,Q exp (iE'x'° - iEx°) 
2m N\ J_ oc 

x {6(x'° - x°) J2 ex P [*(^ - E n)( x '° - x )] W(0, xi) • • • V(0, >4)|n> 

n 

x(n|^(0, Xjv ) • ■ -?(0,xi)|0) + (-1)^0(1° - -T'°)^exp [i(E - E n )(x° - x'°)} 

n 

x(0|^(0,x w ) • ■•^(0,x 1 )|n)(n|V(0,x' 1 ) • • • V(0,x^)|0)} . (2.21) 



Assuming, for simplicity, Eq — (it corresponds to choosing the vacuum energy as the origin of reference) and taking 
into account that 



/oc 
dx°dx'° 8{x'° - x°) exp [-iE n (x'° - a; )] exp [i{E'x'° - Ex ) 
-OC 

= 2irS(E'-E)- 



E-E n + iO 



/OC 
dx°dx'° 6(x° - x'°) exp [-iE n (x° - x m )} exp [i(E'x'° - Ex )} 
-OO 

= -2ttS(E' - E)- 



E + E n - iO 



we obtain 



<?(£) = J2 



$n$r, 



£ - £„ + iO 



is-t-^E 



E + E n - iO 



where the variables x^, ...,x^,Xi, ...,xjv are implicit and 



$ n (x 1 ,...x A ,) = 

S„(x 1 ,...X A r) 



W! 



/ M 



(0|V(0,x 1 ).»^(0,x JV )|n), 
(n|^(0,x 1 )-..V(0,x^)|0>. 



(2.22) 

(2.23) 
(2.24) 

(2.25) 
(2.26) 



In equation (2.24) the summation runs over all bound and continuum states of the system of the interacting fields. 
Let us introduce the functions 



A(.E;xi,...,x^;xi,...,Xjv) = ^<J(#- .E n )$ n (x' 1 ,...,x^ v -)$ n (xi, ...,x N ) 

n 

B(E;x[, ...,x' N ;x 1} ...,x N ) = y^ S(E - E n )S n (x / 1 , ■■■ i x / ^)5 n (xi, ...,x N ) 

n 

These functions satisfy the conditions 

,oo l 

J dEA(E;* 1) ...,x' N] x 1 ,... ) x N ) = —(0\ij(0y i )---ij(0y N ) 



(2.27) 
(2.28) 



f 



dEB(.E;xi,...,x^;x 1 ,...,x A r) = j~ (0\tp{0,x N ) ■ ■ ■ ^(0,xi 



In terms of these functions, the equation (2.24) is 

A(E') 



/"OO 

G{E)= / dE' 
Jo 



x^(0,Xtf)-.-^(0,xi)|0), 

i(0|^(0,x^)---^(0,xi) 

x^0,x' 1 ).^(0,x^)|0). 
B(E>) 



- (-1) N / dE' 

E-E' + iO y ' J E + E'-iO 



(2.29) 
(2.30) 

(2.31) 



where we have omitted the variables Xi, ...,xn,x.[, ...,'x! n and have taken into account that A(E') = B(E') = for 
E' < since -E„ > 0. In fact, due to charge conserv ation , only states with an electric charge of eN contribute to A 
in the sum over n in the right-hand side of equation (2.27) and only states with an electric charge of — eiV contribute 
to B in the sum over n in the right-hand side of equation (2.28). This can easily be shown by using the following 
commutation relations 



[Q,ip(x)] = -eip(x) ) [Q,ip(x)]=eip(x), 



(2.32) 



where Q is the charge operator in the Heisenberg representation. Therefore, the equation (2.31) can be written as 

B(E>) 



Q{E) 



50 dE > A ^ - (-D- r 

(+) E-E' + ii) K ' ./*.<-) 



dE' 

, , E + E 1 -iO 



(2.33) 



where Ei^-1 is the minimal energy of states with electric charge eN and E^,-l is the minimal energy of states with 

mill o./ o rmn o«/ 



electric charge —eN. So far we considered Q{E) for real E. The equation (2.33) shows that the Green function G(E) 

/•OC />00 

is the sum of Cauchy-type integrals. Using the fact that the integrals / dE A{E) and / dE B(E) converge 

J E i+) Je { ~ ) 



(see equations (2.29), (2.30)), one can show with the help of standard mathematical methods that the equation 

Q(E) = r dE 1 4^7 - ("1) W r dE' ^1 (2.34) 

defines an analytical function of E in the complex E plane with the cuts (— oo, E^£\ and [E^^oo) (see Fig. 1). T his 
equation provides the analytical continuation of the Green function to the complex E plane. According to ( 2.33J ) , to 



get the Green function for real E we have to approach the right-hand cut from the upper half-plane and the left-hand 
cut from the lower half-plane. 

In what follows we will be interested in bound states of the system. According to equations (2.24)-( 2~34| ), the 



bound states correspond to the poles of the function Q{E) on the right-hand real semiaxis. If the interaction between 
the electron-positron field and the electromagnetic field is switched off, the poles corresponding to bound states are 
isolated (see Fig. 2). Switching on the interaction between the fields transforms the isolated poles into branch points. 
This is caused by the fact that due to zero photon mass the bound states are no longer isolated points of the spectrum. 
Disregarding the instability of excited states, the singularities of the Green function Q(E) are shown in Fig. 3. The 
poles corresponding to the bound states lie on the upper boundary of the cut starting from the pole corresponding to 
the ground state. It is natural to assume that G(E) can be continued analytically under the cut, to the second sheet 
of the Riemann surface. As a result of this continuation the singularities of Q(E) can be turned down as displayed in 
Fig. 4. 

In fact due to instability of excited states the energies of these states have small imaginary components and, 
therefore, the related poles lie slightly below the right-hand real semiaxis (Fig. 5). However, in calculations of the 
energy levels and the transition and scattering amplitudes of non-resonance processes we will neglect the instability of 
the excited states and, therefore, will assume that the poles lie on the real axis. The imaginary parts of the energies 
will be taken into account when we will consider the resonance scattering processes. 

To formulate the perturbation theory for calculations of the energy levels and the transition and scattering ampli- 
tudes we will need to isolate the poles corresponding to the bound states from the related cuts. It can be done by 
introducing a non-zero photon mass // which is generally assumed to be larger than the energy shift (or the energy 
splitting) of the level (levels) under consideration and much smaller than the distance to other levels. The singularities 
of G{E) with non-zero photon mass, including one- and two-photon spectra, are shown in Fig. 6. As one can see from 
this figure, introducing the photon mass makes the poles corresponding to the bound states to be isolated. 

In every finite order of perturbation theory the singularities of the Green function G(E) in the complex E-plane are 
defined by the unperturbed Hamiltonian. In quantum mechanics this fact easily follows from the expansion of the 
Green function (E — H)^ 1 = (E — H — SV)^ 1 in powers of the perturbation potential SV 

00 
(E - H)- 1 = J2(E - H y l [5V{E - ff rT • (2-35) 

n=0 

As one can see from this equation, to n-th order of perturbation theory the Green function has poles of all orders till 
n+1 at the unperturbed positions of the bound state energies. This fact remains also valid in quantum electrodynamics 
for G{E) defined above. It can easily be checked for every specific diagram in first and second order in a. A general 
proof for an arbitrary diagram is given in Appendix B. 

C. Energy shift of a single level 

In this section we are interested in the energy shift AE a — E a — E a of a single isolated level a of an A^-electron 
atom due to the perturbative interaction. The unperturbed energy Ea is equal to the sum of the one-electron 
Dirac-Coulomb energies 

E^ =e ai +--- + e aN1 (2.36) 



which are defined by the Dirac equation (2.1). In the simplest case the unperturbed wave function u a {x-i, ...,X;v) is a 
one-determinant function 
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U (Xi,...,Xjv) 



W 






{-l) P 4>P ai (Xi) • • ■ V'Pqm (Xiv) , 



(2.37) 



where ip n are the one-electron Dirac wave functions defined by equation (2.1) and P is the permutation operator. In 
the general case the unperturbed wave function is a linear combination of the one-determinant functions 



u a (xi,...,xjv) 



^ C b n -^= V("l) P ^p bl (xi) • • • Vp bN (x w ) • 



b 



'N\ 






(2.38) 



We introduce the Green function g aa (E) by 

gaa(E) = (u a \g(E) 7 1 ... 7 ° N \u a ) 

= / ebq • ■ • dx N dx[ ■ ■ ■ dx' N ul(x[, ...,x' N ) 

x^(-B,x' 1 ,...,x' Ar ;xi,...,x A r)75 ) ---7^u a (xi,...,x A r) 



^From the spectral representation for Q{E) (see equations ( 2.24 )-( 2.34 )) we have 

A, 



9aa(E) 



E-Ea 



terms that are regular at E ~ E a , 



(2.39) 



(2.40) 



where 



Aa = 7VT / dXl ' ' ' dxNdx 'i ' ' ' dx N u l( x i> • ■ ■ > x A')( IV'(0,x' 1 ) • • • ip(0,x' N )\a) 
x(a|V t (0,x JV )..^t(o,x 1 )|0) Ma (x 1 ,. 



•,XjvJ 



(2.41) 



We assume here that a non-zero photon mass /x is introduced to isolate the pole corresponding to the bound state a 
from the related cut. We consider that the photon mass is larger than the energy shift under consideration and much 
smaller than the distance to other levels. To generate the perturbation series for E a it is convenient to use a contour 
integral formalism developed first in operator theory by Szokefalvi-Nagy and Kato [p5H70f . Choosing a contour T in 
the complex E plane in a way that it surrounds the pole corresponding to the level a and keeps outside all other 
singularities (see Fig. 7), we have 



-±-.ldEEg aa {E)=E a A a , 

-^ <f dE g aa {E) = A a . 
2m J r 



(2.42) 



(2.43) 



Here we have assumed that the contour T is oriented anticlockwise. Dividing the equation (2.42) by (2.43), we obtain 

— / dE Eg aa (E) 
2m J T 



E„ 



^-fdE g aa (E) 



(2.44) 



It is convenient to transform the equation ( 2.44 ) to a form that directly yields the energy shift AE a — E a — E a . In 
zeroth order, substituting the operators 



V> in (0,x) = ^ Klpn(x) + ^2 d n^«( X ) ' 



£n>0 



£n<0 



^in(0,x) = J2 6 «^«( x ) + J2 d «^„(x) 



(2.45) 
(2.46) 



e„>0 



e„<0 



into equat ions ( 2.25| ) and (2.26) instead of ip(0,x) and ^(0,x), resp ectiv el y, and cons iderin g the states |n) in (2.25) 
and (2.26) as unperturbed states in the Fock space, from equations (2.24)-(2.26) and (2.39) we find 
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0<°> = 

zaa 



E - 



E (o) 



(2.47) 



JO) 



This e quat ion can also be derived using the Feynman rules for G (see subsection II(E1)). Denoting Ag aa — g aa — j/,,,, 
from ( 2.44 ) we obtain the desired formula |29| 



^- j> dE (E - E^)Ag aa {E) 

1 + ^- <f dEAg aa (E) 
2m ,/r 



AE n 



The Green function Ag aa (E) is constructed by perturbation theory 

Ag aa (E) = Agtt(E) + Ag£(E) + ---, 
where the superscript denotes the order in a. If we represent the energy shift as a series in a 

AE a = A£« + AE& + . . . , 
the formula ( [2.48] ) yields 
1 



(2.48) 



A£« = 



2m 



dEAEAg^(E), 



^-^J^^^(E)-[±J r dEAEAg^E)) [± £ dE A, 



(E) 



(2.49) 
(2.50) 

(2.51) 
(2.52) 



(0) 



where AE = E - E£ 

Deriving equations (2.44) and (2.48) we have assumed that a non-zero photon mass fi is introduced. This allows 
taking all the cuts outside the contour T as well as regularizing the infrared singularities of individual diagrams. In 
the Feynman gauge, the photon propagator with non-zero photon mass (i is 



D pa .(u>,x-y) = -g pa 



or, after integration, 



D pa .(oj,x-y) 



9 P o- 



pa J 


dk 
(2tt) 3 


exp 


(ik- 


(x- 


-y)) 


W 2 - 


-k 2 - 


-M 2 

|x- 


+ i0 


exp (i^/u 2 - 


-M 2 


+ i0 


y|) 



47T X- 



(2.53) 



(2.54) 



where Imy^oj 2 — /i 2 + iO > 0. D pa (ui, x — y) is an analytical function of uj in the complex to plane with cuts beginning 
at the points w = — /U + iO and w = /i — iO (Fig. 8). The related expressions for the photon propagator with non-zero 
photon mass in other covariant gauges are presented, e.g., in |28|. 

As was noted in the previous subsection, the singularities of the two-time Green function in the complex E plane 
arc defined by the unperturbed Hamiltonian if it is constructed by perturbation theory. In particular, it means that 
in n-th order of perturbation theory g aa (E) has poles of all orders till n + 1 at the position of the unperturbed energy 
level under consideration. Therefore, in calculations by perturbation theory it is sufficient to consider the photon mass 
as a very small parameter which provides a separation of the pole from the related cut. At the end of the calculations 
after taking into account a whole gauge invariant set of Feynman diagrams we can put [i — > 0. The possibility of 
taking the limit fi — ► follows, in particular, from the fact that the contour T can be shrunk continuosly to the point 

E = E ( a ] (see Fig. 7). 

Generally speaking, the energy shift of an excited level derived by formula (2.48) contains an imaginary component 
which is caused by its instability. This component defines the width of the spectral line in the Lorentz approximation 
(see sections III(F,G) for details). 

For practical calculations it is convenient to express the Green function g aa {E) in terms of the Fourier transform 
of the 2N— time Green function defined by equation (2.12). By using the identity 



1 

2^ 



dx exp (iojx) = 5(to) 



(2.55) 



12 



one easily finds (see Appendix C) 



where 



g aa (E)6(E - E') 



dp\ ■ ■ ■ dp" N dp'? ---dp 



•$ 



2tt 1 

xS ( E - p° p° N )S(E' -p[° p%) 

x (u a \G{pt ■ ■ ■ .pftrf, • ■ ■ ,P%hi ■ • -7&K) , 



(2.56) 



(u a \G(p' 1 Q ,...,p%;p 1 ,... 1 p N h Q l ... 1 ° N \u a ) 

= / dxi ■ • • dx N dx.[ ■ ■ ■ dx' N u a (x.[, ...,x' N ) 

xG((pf, xi), . . . , (p'l x! N ); (p°, Xl ), • • • , (P%, xjv)) 

X 7l ••■7AT U a(xi,...,X A r). 



(2.57) 



According to equation (2.38) the calculation of the matrix elements in (2.56) is reduced to the calculation of the 
matrix elements between the one-determinant wave functions 



jy. p 



(2.58) 



Uk 



7= 5Z( _1 ) F ^ Ffci ( x ' ' ' V>pfc N ( x jv) 



(2.59) 



To simplify the summation procedure over the permutations in (2.56) which arise from the wave functions as well as 
from the Green function G(p'®, . . . ,p%;p1, ■ ■ ■ ,P%), it is convenient to transform equation ( [2.5(: ) in the following way. 
Denoting G = 0/° ... 7^, we can write 



= X(-i) p G((p'A ! ^i),...,(^V,ep J v);(Pi,ei),-.-,(^,^)), 



(2.60) 



where £ = (x, a) and a is the bispinor index (a — 1,2,3,4). Substituting (2.6C) in (2.56) and using the symmetry of 
G with respect to the electron coordinates, for gtk{E) = (ui\Q(E)-fi ■ ■ ■ 7^1^.) one can obtain (see Appendix D) 



n />oo 

g lk (E)5{E-E>) = -Y J {-l) P r P iM)--^'U N {£,'N) dp\...dpUp'?---dp 

1 p J -00 

x8{E-p\ p%)6(E'-pf P %) 

xG(K°, £), . . . , (p%, &); (p°, fr), • ■ • , (p° n ,Zn)) 

xV'fciKl)---V'feiv(Cjv), 



(2.61) 



where repeated variables {£} imply integration (the integration over x and the summation over a ) . In practical 
calculations by perturbation theory, the formula ( 2.61 ) must be only employed for symmetric sets of Feynman diagrams 
since the symmetry property was used in its derivation. 



D. Perturbation theory for degenerate and quasidegenerate levels 

In this section we are interested in the atomic levels with energies E\ 1 ...,E S arising from unperturbed degenerate 
or quasidegenerate levels with energies E[ , ..., E s . As usual, we assume that the energy shifts of the levels under 
consideration or their splitting caused by the interaction are much smaller than the distance to other levels. The 
unperturbed eigenstates form an s-dimensional subspace $7. We denote the projector on fi by 



P (o) = E p(o) = £ 



UkU' k . 



(2.62) 



fc=i 



fc=i 
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where {uk}% = i are the u npert urbed wave functions which, in a general case, are linear combinations of one-determinant 
functions (see equation ( 2.38J ) ). We project the Green function Q{E) on the subspace fi 



1 (E)^P^g(E) 1 1 ... 1 ° N P^ 1 (2.63) 



where, as in (|2.39), the integration over the electron coordinates is implicit. As in the case of a single level, to isolate 



the poles of g(E) corresponding to the bound states under consideration, we introduce a non-zero photon mass //. We 
assume that the photon mass \x is larger than the energy distance between the levels under consideration and much 
smaller than the distance to other levels. In this case we can choose a contour T in the complex E plane in a way that 
it surrounds all the poles corresponding to the considered states (Ei, ...E s ) and keeps outside all other singularities, 
including the cuts starting from the lower- lying bound states (see Fig. 9). In ad d ition , if we neglect the instability of 



the states under consideration, the spectral representation (see equations ( J2.24|) -( p.34| )) gives 

E-E k 



s t 

g(E) = 2^ -^ j| — I" terms that are regular inside of T , (2.64) 



Using equation (2.64), we obtain 



fc=i 
where 

p fc = P<°>* fc , 4 = $t P (o). (2 . 65) 

As in the case of a single level, in zeroth approximation one easily finds 

p(o) 

5 (0) (^) = Ett^)- (2 ' 66) 

fe=i a - a k 
We introduce the operators K and P by 

K =J_I dE Eg{E) , (2.67) 

2m J T 

P = J_/ dE g {E) , (2.68) 

2m J r 

s 

K = ^Ei<pi<pl, (2.69) 

S 

P = Y,<Pi<fl (2-70) 

t=i 

We note here that, generally speaking, the operator P is not a projector (in particular, P 2 ^ P). If the perturbation 
goes to zero, the vectors {<^i}| =1 approach the correct linearly independent combinations of the vectors {ui\ s i=1 . 
Therefore, it is natural to assume that the vectors {<^i}| =1 are also linearly independent. It follows that one can find 
such vectors {«j}f =1 that 

<p\vk = S ik . (2-71) 

Indeed, let 

s s 

Pi = ^2 at i u 3 ' Vk = 5Z XklUl ■ (2.72) 

3=1 1=1 



The biorthogonality condition (2.71) gives 



/, ciijZkj = Sik ■ (2.73) 

3=1 
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Since the determinant of the matrix {ay} is nonvanishi ng du e to the linear independence of {<fii}f =1 , the system (2.73) 
has a unique solution for any fixed k = 1, ..., s. ^From ( 2.69 )-( [2.71 ) we have 



8=1 
S 

Kv k = ^2 E i<Pi<Pi v k = E k ip k . 



i=l 



Hence we obtain the equation for v k , E k [£9| 

Kv k = E k Pv k . 
According to ( 2.7l| ) the vectors v k are normalized by the condition 

v],,Pv k = 4'fc- 



The solvability of equation (2.76) yields an equation for the atomic energy levels 

dct (K - EP) = . 



(2.74) 
(2.75) 

(2.76) 
(2.77) 
(2.78) 



The generalized eigenvalue problem (2.76) with the normalization condition (2.77) can be transformed by the substi- 
tution ip k = P?v k to the ordinary eigenvalue problem (" Schrodinger-like equation") | p3| 

Hfa = E k ^p k (2.79) 

with the ordinary normalization condition 

1>lil>k> = S kk > , (2.80) 



where H = P"{K)P~^. 

The energy levels are determined from the equation 



det(H-E) = 0. 



(2.81) 



Generally speaking, the energies determined by this equation contain imaginary components which are due to the 
instability of excited states. In the case when the imaginary components are much smaller than the energy distance 
between the levels (or the levels have different quantum numbers), they define the widths of the spectral lines in 
the Lorentz approximation. In the opposite case, when the imaginary components are comparable with the energy 
distance between the levels which have the same quantum numbers, the spectral line shape depends on the process of 
the formation of the states under consideration even in the resonance approximation (see sections III(F,G) for details). 
In what follo ws, ca lc ulatin g the energy levels we neglect the instability of excited states and assume H = (H + H')/2 
in equations ( 2.79 ), ( 2.81 ). 



The operators K and P are constructed by the formulas (2.67) and (2.68) using perturbation theory 

K = K {0) + K (1) + K (2) + ■■■ , (2.82) 

p = p(0) + p(i) + p(2) + ... ( (2.83) 

where the superscript denotes the order in a. The operator H is 

H = H (a) + H (1) + H (2) + ■ ■ ■ , (2.84) 

where 

F<°) = K^ , (2.85) 

(2.86) 



#(1) = #(1) _ lp(l)if(0) _ ijf(0)p(l) 

2 2 



H (2) _ K (2) _ Ip(2)^(0) _ 1 K (0)p(2) 

2 2 

_Ip(!)X (1) - -K {1) P (1) 
2 2 

+-p( 1 )p( 1 )jf (°) + -K {0) P {1) P (1) 



-_p(l)i^(0)p(l) 

4 



(2.87) 
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It is evident that in zeroth order 



k\:> 


=E r 


} 5 lk 


p(0) 


= Sik 




^ 


= E (o 


Sik 



(2.88) 
(2.89) 
(2.90) 



To derive the equations (2.76)-(p.79) we have introduced a non-zero photon mass \i which was assumed to be larger 



than the energy distance between the levels under consideration and much smaller than the distance to other levels. 
At the end of the calculations after taking into account a whole gauge invariant set of Feynman diagrams, we can put 
jU — > 0. The possibility of taking this limit in the case of quasidegenerate states follows from the fact that the cuts 
can be drawn to the related poles by a deformation of the contour T as shown in Fig. 10. 

As in the case of a single level, for practical calculations we express the Green function g(E) in terms of the Fourier 
transform of the 27V-time Green function 

g(E)S(E - E') = ——J dp 1 ...dp N dp' 1 °...dp% 8{E-p\ ■ ■ ■ - p° N )S(E' - pf ■ ■ ■ - p%) 

xP^G(pf, ...p£;p°, ../ Ar ) 7l °...7^ (0) , (2.91) 

where G(pf , —,p%',Pi, ••• ! p°v) is defined by equation ( [2.12; ). 

E. Practical calculations 

In this section we consider some practical applications of the method in the lowest orders of perturbation theory. 
In what follows we will use the notation 

I(uj)^e 2 a p a' T D pa (Lu), (2.92) 

where a p = j a ^ p — (1, a). In addition we will employ the following symmetry properties of the photon propagator 

7(w) = I(-w) , I'(uj) = -I'(-lu), (2.93) 

which, in particular, are valid in the Feynman and Coulomb gauges. Here I' (to) = dl{uj)/ dui . 

1. Zeroth order approximation 



Let us derive first the formula (2.47) using the Feynman rules for G. According to the equations (2.56)-(2.61) we 
have 



g aa (E)6(E -E') = — J2(- l ) P / dx i ' ' -rfxjvdx; ■ • -dx^ 
1 p •* 

/oo 
dp»---dp%dp'?---dp 
-OO 



/0 

N 
— oo 

xS(E - p° p° N )S(E' -pf p%) 

x ^SM, x' l5 X1 )<%1° - p°) • • • 7^S N (p'l x! N , x N )6(p'° - p%) 
X7?---7^a 1 (xi)---^ ajv (x JV ) (2.94) 



= — dp 1 ---dp N dp?---dp> 

* ./-oo 



JY 



x5(E - p° p° N )S(E' - pf p%) 



X- 



2„:- £ „, + »«-«-"ssT^ s «-*. ("«» 



Integrating over the energies one easily obtains equation (2.47) 
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2. One-electron atom 



Formal expressions for the energy shift in the case of a one-electron atom (or in the case of an atom with one electron 
over closed shells) can be derived by various methods. In particular, the Dyson-Schwinger equation can be employed 
for such a derivation. Therefore, the one-electron system is not the best example to demonstrate the advantages of 
the method under consideration. However, we start with a detailed description of this simple case since it may serve 
as the simplest introduction to the technique. 

Let us consider first a diagram describing the interacti on of a one-electron atom with an external potential <$V(x) 
to first order in SV (x) (Fig. 11). According to equation (2.56) we have 



AgV(E)5(E-E') = | J dtdm*^)^ £ Jr^f 



V'ni( x ')V , n 1 ( z ) 

iO) 



^ *x\r< u ,p, p n i y^ V , » 2 (z)V'„ 2 (x) 
i 27T ^^ i — e„ 2 (1 — zO) 



HE 



|ni)(m| 






|n 2 )(n 2 | 



1 



(E - £„) 



£-e ni (l-»0) ^£?_ £na (l_iO) 
a|<JV|a)5(.E' - £). 



|a)<5(£' - £) 



(2.96) 



Substituting ( 2.96 ) into ( 2.51 ), we obtain 



&E^ = ±- J dE {E - e a )AgV{E) 

1 / (a\SV\a) , ,„ , ( 

2m J r E - e a 



(2.97) 



To first order in a, the QED corrections are defined by the self energy (SE) and vacuum polarization (VP) diagrams 
which are shown in Figs. 12 and 13, respectively. Consider first the SE diagram. We find 



A9aa{E)5{E -E) = — dx dydzdx Va( X l^ 2^ Tv~ 



xeiP ^ 6 (E'-p°-u;)^ ^W» 

i 



E' - e ni (1 - iO) 2tt J_ x 

e : niy) t; {z l ^i.,y- z 

*-? V - s„{l - iO) 



ITT 



dw 



xe 7 — 6(p +u-E)—^- — — — — 7> (x) 



E-e n2 (l-iO) 



dw dydz^ a (y)a"J2 



{E-e a y 2tt 
xD pa (uj,y - z)a^ a (z)5(E' - E) 
where a p = 7°7 p = (l,a). Denoting 



V'n(y)V'„(z) 

E - a» - e„(l - iO) 



(a|E(E)|6) 



2- 



dw E 



(an|e 2 a p a CT £> p(T (w)|n&) 



we get 



Affi 1 a , (E) = 



E - u) - e„(l - iO) ' 
a|S(J5)|a) 



(£ - e a f 



(2.98) 



(2.99) 



(2.100) 



Substituting (2.100) into (2.51), we obtain 

Ai#) = ±fdE{E- e a )Ag£(E) 



1 / , (a\T,(E)\a) , , , Nl v 
27T« / r £ - e a 



(2.101) 
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Here we have taken into account the fact that, for a non-zero photon mass (j,, Ag(E) has isolated poles at E = e a in 
every order of pert urbat ion theory. In the final expression one can put fi — > 0. 

The expression ( [2.101 ) suffers from an ultraviolet divergence and has to be considered together with a counterterm 
diagram (Fig. 14). Taking into account the counterterm results in a replacement 

(o|E(e„)|o) -► (o|E*(e„)|a) - (a|(E(e„) - 7 °5m)|a) . 



(2.102) 



The corresponding calculation for the VP diagram (Fig. 13) yields 



2tt 



Ag£)(E)6(E - E') = — r / dxdy ^„(x) 



2tt E' - e a 2?r J_ 



duj' 



9tt 9ir 

xeY—5(E' + J - E)D p Ju',y.- y) e —5{J) 

i I 

i 1 



2tt 



duj Tr 



E 



V'n(y)V' Tt (y) 

w - £„(1 - i0) 



r 



2TrE-e a 



1pa(x) ■ 



Introducing the VP potential by 



we find 



and, therefore, 



tfvp(x) = — / dy a p D pCT (0,x - y) 

27TI ' 



du> Tr 



E 



^(y)^(y) ; 

v - e n (l - iO) 



AffW(£) 



(q|[/yp|a) 



(2.103) 



(2.104) 



(2.105) 



AE ( V = (a\U VP \a) 



(2.106) 

In reality, due to a spherical symmetry of the Coulomb potential of the nucleus, only zeroth components of the a 
matrices contribute to E/vp, 



where 



/i r°° 
dy-. ,/ dwTr[Gc(u>,y,y)) 

^ w-e„(l-»0) 



(2.107) 



(2.108) 



is the Dirac-Coulomb Green function. The expression ( |2. 107 ) is ultraviolet divergent. The charge renormalization 
makes it finite. 

Let us now consider the combined (51^-SE corrections described by the Feynman diagrams presented in Fig. 15. For 
the diagrams "a" and "b" one easily finds 

A<?i 2 a a+6) (£) = fTP \ , 9 ^(a\5V\n)- Er ^ r (n\E(E)\a) 



(E - e a y 

1 
+ {E-e a ) 



E-e n 



•EH s (^)i f 



1 



E-e r , 



(n\SV\a) 



(2.109) 



This contribution is conveniently divided into two parts: irreducible (s n =^ e a ) and reducible (e n = e a ). For the 
irreducible part one obtains 



(n^a) 



J- £ 4E (E - „)*£«**■ \E) = JL £ d E ^L- £ [ 

(a\T,(E)\n){n\6V\a 



(a\SV\n)(n\i:(E)\a) 
E-e n 



E-e n 



(n^a) 



[a\5V\n)(n\Z(e a )\a) 



o|E(e )|n)(n|5V|o) 



(2.110) 
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Here we have taken into account that, due to the spherical symmetry of the Coulomb potential, matrix elements 
(a|E(e Q )|&) are equal to zero if e a = s and a ^ b. The reducible part is 



1 

2ri 



dE (E - e a )Ag, 



(2,a+6,rcd) 



(E) 



1 

2tt7 



dE 



1 



(E - E a ) 



(a\6V\a)(a\Z(E)\a) 



+{a\T,(E)\a)(a\SV\a 

= 2(a\SV\a)(a\X'{e a )\a), 



(2.111) 



where E'(e a ) — (dT,(E) / 'dE) E=e a ■ The re ducib le contribution should be considered together with the related contri- 
bution from the second term in equation (2.52). Taking into account that 



and 



we obtain 



±.£dE(E- e„) Aj,tw»') (E) = (o|«V|a) 

2m J T 2m Jp 



(E-e.y- 



For the total contribution of the diagrams "a" and "b" we find 



a|S'(e a )|o) 



(a|W|a><a|£'(e )|a> 



A£; (2,a+6) _ 



£[■ 



o|«V|n)(n|E(e )|o) , (a\Y>{e a )\n){n\5V\a) 



Sn 



+(a\8V\a)(a\V'(e a )\a) 
For the vertex contribution (the diagram " c" ) we obtain 

1 



&9tt C) (E) 



(E-e a 

><E 



-e 2 — 

) 2 2tt 



du dxdydzi(jl(y)a p ^2——?±— 



An (y)V4 (' 



E-u>-e ni (l-iO) 



SV(x) 



^ E - lu - e n2 (l - i0) 
This diagram is irreducible. A simple evaluation yields 



A 2 (x)V4 2 ( z ) n , x CT , , v 
= —D p<7 {uj,y - zja^ipaiz) ■ 



AE (2,c) 



2tt 



dw / dxdydz ipl(y)a p G c {e a ~ u;,y,x)5V(x) 



xG c (£ a -w,x,z)D po .(w,y- z)a CT <i/> a (z) . 



(2.112) 



(2.113) 



(2.114) 



(2.115) 



(2.116) 



(2.117) 



The rel ated m ass-counterterm diagrams (Fig. 16) are accounted for by the replace ment S — > E# = E — 7°<5m in 
equation (2.115). This replacement makes the irreducible contribution in equation (2.115) to be finite. As to the 
reducible contribution, its ultraviolet a nd inf rared divergences are cancelled by the corresponding divergences of the 
vertex contribution given by equation (2.117). 



3. Atom with one electron over closed shells 

The consideration given above can easily be adopted to the case of an atom with one electron over closed shells by 
regarding the closed shells as belonging to a redefined vacuum. The redefinition of the vacuum results in replacing 
iO by — iO in the electron propagator denominators corresponding to the closed shells. In other words, it means 
replacing the standard Feynman contour of integration over the electron energy C by a new contour C" (Fig. 17). 
In this formalism the one-electron radiative corrections are incorporated together with the interelectronic-interaction 
corrections and the total energy of the closed shells is considered as the origin of reference. The difference of the 
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integrals along C and C is an integral along the contour C- m t . It describes the interaction of the valence electron with 
the closed-shell electrons. Therefore, to find the interelectronic-interaction corrections we have to replace the contour 
C in the expressions for the one-electron radiative corrections by the contour C- ln f For example, in the case of one 
electron over the (Is) 2 shell in a lithiumlike ion, the first order interelectronic-interaction corrections are obtained 
from the formulas for the SE and VP corrections derived above by the replacement 

n v ' c 

As a result of this replacement, one obtains for the interelectronic-interaction correction 

(ee=ei.) 

AjB (i,int) = J2 [{ac\I(0)\ac)-{ac\I(e a -e l8 )\ca)}, (2.119) 



where I(u>) is defined by equation (2.92). In Ref. |3q| this formalism was employed to derive formal expressions for 



the interelectronic-interaction corrections to the hyperfine splitting in lithiumlike ions (see section IV(C2)). 

4- Two-electron atom 

Let us consider now the energy shift of a single level (n) in a two-electron atom. To first order in a, in addition to 
the one-electron SE and VP contributions (Figs. 18,19) we have to consider the one-photon exchange diagram (Fig. 
20) . The derivation of the energy shift due to the SE and VP diagrams is easily reduced to the case of a one-electron 
atom by a simple integration over the energy variable of a disconnected electron propagator. Therefore, below we 
discuss only the one-photon exchange diagram. 

For simplicity, we assume that the unperturbed wave function of the state under consideration is the one-determinant 
function 

u„(xi,x 2 ) = -^=^(-l) P 7/> Pa (xi)?Ap b (x 2 ). (2.120) 



The transition to the general case of a many-determinant function (2.38) causes no problem and can be done in the 
final expression for the energy shift. 



According to equation (2.56), for the one-photon exchange diagram we have 

2 r°° ^^ 1 1 



• n /*oo 



p'O _ =.„ _i_ „-n p _ ^'0 
1 1 



; , PT -ep a +iQ E- pf - e Pb + z0 



' Pi - £ a + iO E - p\ - E b + i0 



(PaPb\I(pT - p")\ab) . (2.121) 



Formula (2.51) gives 

i r ' \ 2 poo 



r ~ \2ttJ J^' l ^ l ^ p'f-epa+iO 

1 1 1 



E - pf - e Pb + iO pi - e a + iO E-p\-e b + iQ 
x{PaPb\I($f-p\)\ab), (2.122) 



where, as in (2.51), AE = E — En ' . Transforming 



(0) 



1 1 1/1 1 



u. 



p[ Q -epa + iOE- p[° - ep b + iO ~ AE \pf -£ Pa + iO E- p[° - s Pb + iO 

1 1 _ 1 / 1 1 

p\ - e a + iO E - p\ - e b + iO ~~ ~AE \p\ - e a + iO + E - p\ - e b + iO 

we obtain 
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(2.123) 

(2.124) 






(2.125) 



The expression in the curly braces of ( p. 125 ) is a regular function of E inside the contour T, if the photon mass \x 
is chosen as indicated above. A direct way to check this fact co nsist s in integrating over p\ and p'® by using the 
apparent expression for the photon propagator given by equation ( 2.53| ). It can also be understood by observing that 
the integrand in this expression is the sum of terms which contain singularities from the electron propagators in p® 
(p'l) only above or only below the real axis (for real E). Therefore, in each term we can vary E in the complex 
E plane within the contour T, keeping the same order of bypassing the singularities in the pi (p'i) integration by 
moving slightly the contour of the p\ (p'f) integration into the complex plane. (In contrast to that, the contour 



of the pi (p'i) integration in equation (2.122) is sandwiched between two poles and, therefore, can not be moved 
to the complex plane.) The branch points of the photon propagators are moved outside the contour T due to the 
non-zero photon mass. However, instead of investigating the analytical properties for every specific diagram, it is 
more convenient to use the general analytical properties of the Green function g nn (E) in the complex E plane in 
every order of perturbation theory (see the related discussion above and Appendix B). According to these properties, 

the function Ag n7 l(E) can no t have poles at the point E = En of order higher than 2 and, therefore, the expression 



in the curly braces of ( |2.125 ) is a regular function of E at this point. So, we have to calculate the first order residue 
at the point E = E n ' . We also stress that we do not need any apparent form for the analytical continuation of the 
expression in the curly braces to the complex E plane since we calculate it only for real E, at the point E 
where the present expression is valid. Calculating the E residue we obtain 



p(o) 



■ 2 /*oo 

A£; « 1) = (f) / *«£(-i 

V»9-e„-H0 



1 



1 



P? 



(p? - e a ) + tO 



Taking into account the identity 



\x + 



1 



2tt \x + iO —x + iO 



-ep a + iO -(Pi-spa) 
(PaPb\I(p[ -p1)\ab). 

= 6{x) 



iO 



we find 



AM 1 



p 



(-l) p (PaPb\I(sp a -e a )\ab) 



(2.126) 



(2.127) 



(2.128) 



5. General rules for practical calculations 

Let us consider now some general remarks to the derivation given above. In order to perform first the integration 
over E we have separated the singularity in AE by emp loying the i dentitie s ( [2.125 ) and ( p. 124 ) . Another way could 
consist in transforming the left-hand sides of equations (2.123) and (2.124) by the identity 



2tt 



p° - s a + »0 E - p° - s b + iO 



S(P° 



1 



p° - e a - iO E - p° - e b + iO ' 



(2.129) 



where we have used equation ( 2.127 ). Using this identity allows one to separate contributions singular in 1/A E from 
non-singular ones. The singular contributions result only from the first term in the right-hand side of equation ( |2.129 ). 
It can easily be understood by observing that the second term has both singularities in p° above the real axis (for 
real E). It follows that the contour of the p° integration in the expression for the energy shift can be moved slightly 
into the complex E plane keeping the same order of bypassing the singularities. It means that we can vary E in the 
complex E pla ne wit hin the contour T and, therefore, the integrand is a regular function of E within this contour. 
Identities like (2.129) are very useful in calculations of three- and more electron atoms (see section IV(B3)). 
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We also want to note that in all cases the order of the singularity in 1/AE is quite evident from the type of the 
diagram under consideration. If the diagram is irreducible, the factors 1/AE ma y com e only from the initial and 
final propagators. In this case the second term in the right-hand side of equation (2.129) does not contribute to the 
energy shift and, therefore, the derivation of the formal expression for the energy shift becomes trivial. For reducible 
diagrams the factors 1/AE arise also from internal electron propagators. 

W e can for mulat e th e follo wing simple rule for deriving the energy shift from a certain diagram. Using identities 
like ( [2.123 ), (2.124) or (2.129), we separate all singularities in 1/AE and then integrate over E assuming that the rest 
is a regular function of E within the contour T . As is discussed above, the order of the singularity is quite evident 
for every specific diagram and it is a simple task to separate the factor 1/AE to the right power. However, even if 
we separate this factor to a power which is larger or smaller than the real order of the singularity, it is impossible to 
miss the correct result. In the first case (the power is larger than the real order of the singularity) the result of the 
calculation remains the same as in the case when we separate the factor 1/AE to the right power. In the second case 
(the power is smaller than the real order of the singularity) we obtain an infinite result (~ 1/0). It means that we 
should increase the power of the separated singularity and repeat the calculation until we get a finite result. 



6. Two-photon exchange diagrams for the ground state of a heliumlike atom 



The two-photon exchange diagrams are presented in Fig. 21. Here we derive the energy shift from these diagrams 
for the case of the ground state of a heliumlike atom. The case of a single excited state of a two-electron atom is 
considered in detail in 1 36 1 . The wave function of the ground state is given by 



»i 



(Xi,X 2 ) = ^=^(-l) P V'Pa(xi)'(/'Pfc(x2). 



(2.130) 



The unperturbed energy is E\ = e a + £b, where e a = Eb 



Consider first the two-photon ladder diagram (Fig. 21a). For the first term in (2.52J) we have 



^ ( 2 = ^/^b-i>'(^ 



2n 



dp\dp'°duj 



J2 {PaPb\I( P '° - w)|nina)<nin 2 |/(w - P°)\ab) 



1 



1 



1 



' pf -ep a +iOE- pf -e Pb + iQio- e ni (1 - iO) 
1 1 1 



E-w- e„ 2 (1 - iO) p\ - e a + iO E - p\ - e b + iO ' 
Let us divide this contribution into irreducible (e ni + s n2 ^ e a + e b ) and reducible ( e ni + e n . 3 = e a + Eb) parts 



(2.131) 



.(2) 



A^ - A£ lad 



(2,irrod) 



A£, 



(2, red) 
lad 



(2.132) 



Using the identities (2.123) and (2.124) we obtain for the irreducible part 



JLfdE AEA 9 [f^\E) = ±-J d E -^{E(-D P (i X 3 



27ri 



2-ni 



A£H^ V 



V2tt 



dp 1 dp 1 dio 



^n 1 +£n 2 ¥ zl ^a+eb 

x Y, <PaP6|/(p , 1 -w)|niJi2)(nin 2 |/(cj-p°)|a6) 



711,712 



( L 

'VrfP-ep, 



1 



^p'l>~ep a + iO ' E-pf-epb + iQ 

1 1 

uj - e ni (l-iO)E-u- e n2 (1 - iO) 



U?-e, 



Pi - e a + iO E- p\ -£ b + iO; 



(2.133) 
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The expression in the curly braces of (2.132) is a regular function of E inside the contour T if the photon mass p is 
chosen as indicated above. (If it were not so, we would get an infinite result; see the related discussion in the previous 
subsection.) Calculating the E residue we find 



Ai4 2 f cd) = E(~!) P ^ J ^ E (PaPb\I(e Pa W )| 



x(nina|J(w-e )|o6) \, -r -^ l — -. (2.134) 

w - s ni (1 - «0) E (°) _ ^ _ e„ 2 (1 - iO) 

This derivation shows that the energy shift from an irreducible diagram is obtained by evaluation of the " usual S- 



matrix" element. For the numerical evaluation of (2.134) it is convenient to rotate the contour of the integration in 
the complex to plane |7l|] . 

For the reducible contribution we have 

x J2 (PaPb\I(p' 1 -Lu)\n 1 n 2 )(n 1 n 2 \I(uj-p 1 )\ab) 

ni,ri2 

X \pf - e Pa + !(j + E-p^-epb + io) 

x ( I | 1 ) 

\uj — e ni + iO E — lu — e n2 + iO)/ 

X \p° 1 -e a + iO + E-pi-et + io)}- (2 ' 135) 



The expression in the curly braces of ( 2.13 5| ) is a regular function within the contour T. Calculating the E residue 



and taking into account that e a = £{,, we get 

£ nl +e„ 2 =2s a 



^-fdEAEAg£ Icd) (E) = -^-Y.(- l ) P E {/ dp[° (PaPblKp? - e a )\ ni n 2 



x(mn2|I(£ a -Pi)|a6) 



«1,«2 

1 

( :„ - pf + ioy 
1 



x{mn 2 \I(0)\ab)- , - + / drf (PoP6|J(0)|nina> 



(e a - p? + zO) 2 
+ / dw (PoP6|J(e -w)|nin 2 ) 



{ ni n 2 \I{u-e a )\ab)- ^—^j. (2.136) 

e„ — u> + iv) z j 



This contribution should be considered together with the second term in equation (2.52). As was obtained above, the 
first factor in this term is 

— (f dE AEAg^fE) = V (-if (PaPb\I(0)\ab) . (2.137) 

2m Jr p 

A simple calculation of the second factor yields 

±-J t *E A 9 i!»(E, . ±£4E ^L ? {(^) 2 £ #Wl 

1 1 



V(-D' 

^ V ' \p'l'--e a + i0 • E-p'l>-e a +i0 
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+ J-oo dPl (/- £a -»0)2 (PflPb|/(P °~ £a)|ab) }- (2 ' 138) 



For the total reducible contribution we obtain 



AE ^ ' = ti% dE AE A ^ ' w - Us £ di? Ai? A5 " <*>; fe A di? A5 " w 

e ni +e n2 =2e a . 



1 

(w - e"T- iO) 2 



x(nir»2|/(w-ea)|o6>7^^: 7K^ ■ ( 2 - 139 ) 



A similar calculation of the two-photon crossed-ladder diagram (Fig. 21b) gives 

/oo 
doj (Pan 2 \I(oj - e a )\nib)(nxPb\I(w - e a )\an 2 ) 

*o;-e ni (l-iO)a;-e n2 (l-iO)" (2 ' 140) 

The contribution A_Ej a( j IC contains an infrared divergent term which is cancelled by a related term (e ril = e„ 2 = e a ) 

(2) 

from the contribution AE'cross- In the individual contributions the infrared singularities are regularized by a non-zero 
photon mass /i. If the ladder and crossed-ladder contributions are merged by the common ^-integration, the integral 
is convergent and we can put /i — ► before the integration over lo (see J36j for details). However, to show how the 
calculation for a non-zero photon mass can be performed, let us calculate the reducible contribution for a finite /i. 
We have to calculate the integral 



j_ f°° duj exp [i^io 2 - /i 2 + ^0 (n 2 + r 34 )] 1_ 

2?r 7.^ ri 2 r 34 (lo - iO) 2 



h = J_l duj ^L'V" _^ J'"^ 12T ^ ; ^ . (2.141) 



Using the identity 



we obtain 



exp 



[*^ 2 -M 2 +*0 r]=-- J" dk k . 2 S ,^ (fcr) 2 - , (2.142) 



7T 







' (cj 2 - jfc 2 - m 2 + tO) ' 



7T 2^ r 12 r 34 Loo 7o (^ 2 - fc 2 - M 2 + i0) (w - *0) 2 



^=-Z^Z^— / <W ^ ,..2,2 ..2 -nw.. ^ 2 - (2.143) 



Decomposing the denominator 



lo 



2 1.2 



k 2 - p? + i0 = (lo - Vfc 2 + /1 2 + i0)(w + VF+7? - *0) (2.144) 



and integrating over lo, we find 



,,,-l^r^ siffiraj^l) . (2 , 45) 

7rri2r34y ( fc + m) 7 



According to |72], the last integral is 



where 



/1 = -- ri2+r34 ivo[ M (r 12 + r 34 )] , (2.146) 

7T r i2 r 34 
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K (z) = log (z/2) £ %f- + E M^( fc + 1) • ( 2 - 147 ) 

fc=0 V '' k=0 v '' 



Considering fi — * wc find 

A^ ( a 2 d r ° d) =--E(-!) P /rfx 1 ...dx 4 ^ Pa (x3)?p b (x 4 ) 73 V 

x V ^ ni (X3)^ ni (xi)?/>„ 2 (x 4 )V'„ 2 (x 2 )7i A 72 ( — + — ) 

Vrl2 r34/ 

x[log(n 2 +r 34 ) +log^-log2-^(l)].g pCT g Al/ 7/' a (xi)-0 b (x2). (2.148) 

The corresponding contribution (e ni = e n2 = e a ) from the crossed-ladder diagram is calculated in the same way. 
The sum of the reducible contribution from the ladder diagram and the related contribution from the crossed-ladder 
diagram is [|9[ 

A£ (2 ' infr) = -— £(-l) P E' / dxi ' ' ' dx4 ?Pa(x 3 )?P & (x 4 )73 P 74 
x-0 ni (x3)V> ni (xi)i/>„ 2 (x 4 )i/>„ 2 (X2)7 4 A 72 

x S^S^ ( 1 ) log (ri2 + r 34 ) 

L Vr i2 r 34 / 

-9 pv 9a\( 1 ) log(ri 4 + r 23 ) Va(xi)V>&(x2). (2.149) 

Vri4 r 23 / J 

The terms containing the factor log/i — log 2 — ip(l) have cancelled each other. For the numerical evaluation of this 
contribution, one can use the following representation |3q ] 

AjE (2,xnfr) = I rdy {ba \ S '(y)\ab)[2(ab\S(y)\ab) - (oo|S(y)|oo) - (bb\S(y)\bb)] 

n Jo y 

+ 1 r *Z (ba\S(y)\ab)[2(ab\S'(y)\ab) - (aa\S f (y)\aa) - (bb\S' (y)\bb)} 
k Jo V 

-1 r <*L (ba\S(y)\ab)(ba\S'(y)\ab) , (2.150) 

Wo y 

where 

= a (l-ar a2 ) ryii (2 151) 

ri2 
5'(y) = -a(l-o; 1 -a2)e~ yri2 . (2.152) 

Other representations of this term can be found in Refs. [ |7l|j73|| . 

7. Self-energy and vacuum-polarization screening diagrams 

The self-energy and vacuum-polarization screening diagrams are presented in Figs. 22 and 23, respectively. The 
derivation of the calculation formulas for the energy shift of a single level due to the SE screening diagrams using the 



TTGF method is described in detail in Ref. []43j. If the unperturbed wave function is given by equation (2.12C), the 
contribution of the SE screening diagrams is written as 

AE = A£ irrcd + AE Icd + AE VCI , (2.153) 

A£ irred = E(-l) P { E (PaPb\I(*)\nb)^—{nMe a )\a) 

P n £a Sn 



V (PaPb\I(A)\an) (n|E(e 6 )|6) 
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- Y (Pa\X(e Pa )\r 



£Pa — £n 



(nPb\I(A)\ab) 



En^EP 



T (Pb\Z(e Pb )\n) (Pan\I(A)\ab)\ , 

'-~ l £Ph — Sr, J 



AP rcd = (ba\I'(E b - e a )\ab) [(o|E(e„)|a) - <6|S(e 6 )|6> 

+ ^(-l) p (PaP&|/(A)|a6)[< a |S'(e a )|a) + (b\Z'(e b )\b) 



AE™ = J>1) P Y 



i 
2^ 



dui 



{ ni Pb\I{A)\n 2 b) (Pon 2 |/(w)|nia) 



£ Pa - w - £m (1 - «0)] [e a - u - e„ 2 (1 - iO)] 

(Pan 1 \I(A)\an 2 )(Pbn 2 \I{oj)\n 1 b) 



[e Pb - w - e ni (1 - *0)] [ei, - w - e„ 2 (1 - iO)] 
where A = ep a — e a . The corresponding contribution of the VP screening diagrams is 

AE = AE'^ cd + AE r a cd + AE b , 

AK rred = £(-l) P { Y {PaPb\I{A)\nb)-±— (n\U VP \a) 



- J2 (PaPb\I(A)\c 

n 

Y (Pa\U VP \i 



1 



1 



(n|[/ V p|6) 



nPb\I(A)\ab) 



n ¥=ep b 

V (P6|J7vp|n) (Pan\I(A)\ab) \ , 

z — ' eph — e n J 



£Pa — £• 

(P6|C/ V p|n! 

epb — e 

AP- d = (ba\I'(e b - e a )\ab) [<o|J7vp|o) - (b\U VP \b) 

AE b = ^(-l) p (PaP6|/ V p(A)|a6) , 



(2.154) 



(2.155) 



(2.156) 



(2.157) 



(2.158) 

(2.159) 
(2.160) 



where 



a" I I I ai M exp(i|e||x-z 1 |) a 2 „exp(«|e||y - z 2 | 
Jvp (e, x, y) = — / da; j dz 1 \ dz 2 — ' 



X-Zi 



xTr[a»G(uj-^z l ,z 2 )a v G(u J +^z 2 ,z 1 )] 



|y-z2i 



(2.161) 



The expressions ( 2.153 )-( 2.16C ) are ultraviolet divergent. The renormalization of these expressions can be performed 
in the same way as for the first-order SE and VP contributions (see Refs. p9|,p0|,p2|,p3[ for details). 



8. Quasidegenerate states 

Let us now consider some applications of the method to the case of quasidegenerate states of a heliumlike ion. This 
case arises, for instance, if one is interested in the energies of the (ls2p 1 / 2 )i and (ls2p 3 / 2 )i states. These states are 
strongly mixed for low and middle Z and, therefore, must be treated as quasidegenerate. It means that the off-diagonal 
matrix elements of the energy operator H between these states have to be taken into account. The unperturbed wave 
functions are written as 



Ui(xi,x 2 )=A N Y {3ti m ti3t2 m i2\ JM )Y^ l ^ P ^ Pl ^ l ^ p, ^ yi2 ^ 



rrii, m 



(2.162) 
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where An is the normalization factor equal to l/v2 for non-equivalent electrons and to 1/2 for equivalent electrons, 
J is the total angular momentum, and M is its projection. However, in what follows, to compactify the formulas we 
will construct the matrix elements of H between the one-determinant wave functions 



j(x 1 ,x 2 ) = - ? =^(-l) P V'p ll (xi)V'p 42 (x2) 



(2.163) 



The transition to the wave functions defined by equation ( 2.162J ) can easily be accomplished in the final formulas 



(1) 



First we consider the contribution from the one-photon exchange diagram. To derive the formulas for H ik we will 

assume that E\ ' ^ E k ' . However, all the final formulas remain to be valid also for the case E\ ' = E\' which 
was considered in detail above. According to the Feynman rules and the definition of g(E), the contribution of the 
one-photon exchange diagram to g^'(E) is 



fl « (2S) = (£) _ fcWEc- 1 ) 



i 
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Pi -e kl +iQE-Pi- e k2 



Vi - £Ph +tOE- pf - e Pl2 
{PiiPi 2 \I{Vi-vi)\kik2). 



i() 



£0 



(2.164) 



Using the identities (2.123) and (2.124), we obtain 
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< " 53 * JE 
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iO, 



(2.165) 



The expression in the curly braces of (2.165) is a regular function of E inside the contour r, if the photon mass fi is 



chosen as indicated above. Calculating the E residues and taking into account the identity (2.127), we obtain 
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(2.166) 



In the same way we find 
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Symmetrizing equations (p. 166) and (2.167) with respect to both electrons we transform them to the form 



(2.167) 
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.(w + Ai - iO)(w - A 2 - iO) (w + A 2 - iO)(w - A x - iO). 
where Ai = ep^ — ej-j and A 2 = en 2 — £k 2 . Substituting ( 2.168 ), ( 2.169J ) into ( 2.86J ), we get 

H% = \ ^(-l) P [(PiiPz 2 |/(A 1 )|fc 1 fe 2 ) + (PiiPt 2 |I(A 2 )|fc 1 fc 2 >] . 



(2.168) 



(2.169) 



(2.170) 



Let us now consider the contribution to H from the combined SV - interelectronic interaction diagrams presented in 
Fig. 24. For simplicity, we will assume that SV is a spherically-symmetric potential. In the case under consideration, 
the simplest way to derive the formulas for H ik ' consists in using the fact that these diagrams can be obtained as the 
first-order (in SV) correction to the one-photon exchange contribution deriv ed above. So, the contribution from these 
diagrams can be obtained by the following replacements in equation ( 2.17CJ) 



where, to first order in SV, 



|*i)-»|fci)+*|fci), 
\k 2 ) -» |fo) + S\k 2 ) , 

{Pi!) ^ \PH) + S\Ph) , 
\Pi 2 ) -> |Pi 2 ) + 5\Pi 2 ) , 
7(e a - £&) -» 7(e a + Se a -e b - 5e b ) , 



fc a = (a\8V\a) , 

" a \n)(n\SV\a) 



S\a) = J2 



£a-£ 



(2.171) 
(2.172) 
(2.173) 
(2.174) 
(2.175) 



(2.176) 

(2.177) 



a °n 



Here we have taken into account that, due to the spherical symmetry of SV, (n\SV\a) — if e n — e a and \n) ^ |o) . 
Decomposing the modified expression for the one-photon exchange diagram to the first order in SV, we find that the 
total correction is the sum of the irreducible and reducible parts, 



n ik 
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(2.178) 



where 



H%™ d) = ^(-l) p [(^iPz 2 |/(A 1 ) + I(A 2 )\k 1 k 2 



-(Pi 1 SPi 2 \I(A 1 ) + I(A 2 )\k 1 k 2 ) 
-{Pi 1 Pi 2 \I{A 1 ) + I{A 2 )\Sk 1 k 2 ) 
-(Pi 1 Pz 2 |7(A 1 ) + 7(A 2 )|fc 1 5fe 2 ) 



(2.179) 



and 



#i 2,r ° d) = ^E(- 1 ) P ^<Piil^|Pii) - (fc 1 |<5T/|fe 1 )](Pi 1 P« 2 |/ , (A 1 )|fe 1 fc 2 
p 
+[<P» 2 |W|Pi 2 ) - (fe 2 |^|fe 2 )](Pi 1 Pi 2 |7'(A 2 )|feife 2 )}. 



(2.180) 
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Equations (|2.179| ) and (2.18C) provide the matrix elements between the one-determinant wave f unctio ns defined 
by equation ( |2.163| ). To get the matrix elements between the wave functions defined by equation ( [2.162 ), we have 
to multiply these equations by the Clebsch-Gordan coefficients and to sum over the projections of the one-electron 
angular momenta. 

(2) 

The expression for H ik can also be derived by direct application of the TTGF method. It can easily be performed 

(2) 

in the same way as for the one-photon exchange diagram. We note that this derivation yields a formula for H ik 
which is slightly different from the expression given above. In particular, for the irreducible contribution one finds 
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(2.181) 
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The term AH, 

which can be represented as 
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(2.182) 
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approaches zero if P) — » E k ' . The expressions (2.179) and (2.181) differ by the term AiP 
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Here O' 2 ' is an operator of second order in the perturbation parameter which we denote by Ao (for simplicity, we 
assume here that 6V and the interelectronic-interaction operator I(lo) are characterized by the same perturbation 
parameter). This fact can be understood by observing that the integrand in (2.182) is the sum of terms which contain 
singularities from the external electron propagators only above or below of the real axis and, therefore, integrating over 



lo cannot result in the appearance of contributions ~ 1/(P,- P k ) that could compensate the factor (E\ 



(2) 

particular, it means that O ik remains finite when E, 



M 



E 
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It can be shown that the term AH 



(2,irrcd) 



ED- In 
contributes 



only to third and higher orders in Ao and, therefore, can be omitted if we restrict our calculations to second order 
in Ao. Let us prove this fact for the case of two quasidegenerate levels. In this case the energy levels are determined 
from the equation 



(E - H n )(E - H 22 ) - H X2 H 2l = 



(2.184) 



2<.) 



which yields 



£1,2 = — I g ' 2 ± ^ V(gn - H22) 2 + 4ff 12 g 21 . (2.185) 



If E\ — E2 ~ Ao, the proof of the statement is evident. If E\ — E2 3> Ao, the contribution of the second-order 
off-diagonal matrix elements is given by 



AE l>2 « i^tf^ + H$H%>)/(EW - 25<°>) . (2.186) 

om this equation we find ' 
Ao. 



^From this equation we find that the terms ~ (E\ — Ej_ )O ik in H ik contribute only to third and higher orders in 



Equations (2.17S )-(2.180) will give a part of the VP screening contribution, if we replace SV by t/yp- The total 
contribution of the VP screening diagrams for quasidegenerate states is derived in Ref. p4fl . Corresponding formulas 
for the SE screening diagrams are derived in Ref. [ J75| |. 

F. Nuclear recoil corrections 

So far we considered the nucleus as a source of the external Coulomb field Vq ■ This consideration corresponds to 
the approximation of an infinite nuclear mass. However, high precision calculations of the energy levels in high-Z 
few-electron atoms must include also the nuclear recoil corrections to first order in m/M (M is the nuclear mass) and 
to zeroth order in a (but to all orders in aZ) . As was shown in |4l| , these corrections can be included in calculations 
of the energy levels by adding an additional term to the standard Hamiltonian of the electron-positron field interacting 
with the quantized electromagnetic field and with the Coulomb field of the nucleus Vc- In the Coulomb gauge and 
the Schrodinger representation, this term is given by 

Hm = ^ |dx^( x )(-jV x )^(x) [ dy i>Hy)(-iV y )i>(y) 

~ e -§Jd* ^(*)HV x ty(x)A(0) + ^A 2 (0) . (2.187) 

Hm taken in the interaction representation must be added to the interaction Hamiltonian. It gives the following 
additional lines and vertices to the Feynman rules (we assume that the Coulomb gauge and the Furry picture are 
used). 

1. Coulomb contribution. An additional line ("Coulomb-recoil" line) appears to be 

^%f°° du. 

x y 

This line joins two vertices each of which corresponds to 



W3 



(j 1 




uJ 2 



-27rJ7°(5(wi - u>2 — uj 3 ) J dxpk 
where p = — iV x and k = 1, 2, 3. 
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2. One-transverse-photon contribution. An additional vertex on an electron line appears to be 

-27ri7°<5(wi -oj 2 -uj 3 )^ J dxp k , 
The transverse photon line attached to this vertex (at the point x) is 

x y 

At the point y this line is to be attached to a usual vertex, where we have — 2irie~f ai2ir5(u>i — u)2 — 0J3) J dy. 
oei (I = 1, 2, 3) are the usual Dirac matrices. 



3. Two-transverse-photon contribution. An additional line ("two-transverse-photon- recoil" line) appears to be 

: - -gr I-oo duj D «(w, x)Ak(w, y) ■ 



. . i e 2 Z 2 



x y 

This line joins usual vertices (see the previous item). 

Let us apply this formalism to the case of a single level a in a one-electron atom. To find the Coulomb nuclear 
recoil correction, we have to calculate the contribution of the diagram shown in Fig. 25. According to the Feynman 
rules given above we obtain 

Ag W (E) = * JL± r ok, y (aNn) / nN " } • (2-188) 

AEc = JL A r dcJ y (°W"X"W«) ( lg9) 

The one-transverse-photon nuclear recoil correction corresponds to the diagrams shown in Fig. 26. A similar calcu- 
lation yields 

■iiraZ i f 00 ^ J (a\pi\n){n\a k D ik {e a - u)\a) 



The formula ( 2.51 ) gives 
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(o|afcAfc(eo-w)|n>(n|pi|a> 



w-s n (l -iO) j 

The two-transverse-photon nuclear recoil correction is defined by the diagram shown in Fig. 27. We find 

(47raZ) 2 i 
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The sum of all the contributions is 

1 i f°° 
AE =T7TT duj (°I[p* + 47roZojAi(w)] 

xG c (w + e a )[p, + 4iraZa m D mi (uj)}\a) , (2.192) 
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where Gc(w) is the Dir ac- Cou lomb Green function defined above. For practical calculations it is convenient to 
represent the expression ( 2.192 ) by the sum of a lower-order term A£l and a higher-order term AEr: 



AE = AE L + AE H , 
1 



AE 



l = 2M<«l[Pi " (A(O)pi + Pi D i (0))]\a) , 



AE H = 



2ttM 



dio 



(a|(A(w) 



\Pj,Vc] 
ui + iO 



Gc(u + e a )(pi(u) 



ui + iO 



(2.193) 

(2.194) 

(2.195) 



where D^u) = —AiraZaiDi^oj) . The term AE^ contains all the recoil corrections within the (aZ) A m 2 jM approx- 
imation, while the term AE-r contains the contributi on of or der (a Z) 5 m 2 /M and all contributions of higher orders 



in aZ which are not included in AE-^. The formulas (2.192)-( 2.19q) w ere first derived by a quasipotential method in 
Ref. [p0| and subsequently rederived by other methods in Refs. ]76|,[77[. The representation (2.192) was found in Ref. 



Consider now a two-electr on ato m. For simplicity, as usual, we assume that the unperturbed wave function is 
a one-determinant function ( 2.120| ). The nuclear recoil correc tion is the sum of the one-electron and two-electron 
contributions. Using the Feynman r ules an d the formula ( 2.51 ), one easily finds that the one-electron contribution is 
equal to the sum of the expressions (2.192) for the a and b states. The two-electron contributions correspond to the 
diagrams shown in Figs. 28-30. A simple calculation of these diagrams yields 



A£ (int) = ^^(-l) P ( p «l \pi+iiraZai D u (e Pa -e a 



)}\a) 



x (Pb\\pi + <±iraZa m D mi (sp b - s b )]\b) 
The formula (2.196) was first derived by a quasipotential method in Ref. [plf . 



(2.196) 



III. TRANSITION PROBABILITIES AND CROSS SECTIONS OF SCATTERING PROCESSES 

According to the basic principles of quantum field theory fl7q] , the number of the particles scattered into the interval 
dp[ ■ ■ ■ dp' r for a unit time and in a unit volume is 



dW, 



p^p> 



(2^) 3s+1 m • • • n s \T p , p \ 2 6(E p - E p ,)dp[ ■ ■ ■ dp' r 



(3.1) 



where p, p' are the initial and final states of the system, respectively; t p > p is the amplitude of the process defined by 



{p'\(S-I)\p)=2niS(E p -E pl y 



p p ' 



(3.2) 



S is the scattering operator, / is the identity operator, s is the number of initial particles, r is the number of final 
particles; n%, ...n s are the average numbers of particles per unit volume. 

We will consider the scattering of photons and electrons by an atom that is located at the origin of the coordinate 
frame. The differential cross section is defined by 



da 



dWr, 



J 



(3.3) 



where j is the current of initial particles (for photons j — nc; for electrons j — nv, where v is the velocity of the 
electrons in the nucleus frame). The total cross section can be found by integrating the differential cross section over 
all final states. The cross section for the elastic scattering is 



(clast) 



(P) 



(27T) 



■ j dp' 8{E p -E p ,)\t p , p \ 



The total (elastic plus inelastic) cross section can be found by employing the optical theorem 

CTtot(p) = IniTpp . 



(3.4) 



(3.5) 



In terms of the amplitude f p / p which is defined so that da — \f p ' p \ 2 dfl, the optical theorem has a well-known form 
(see, e.g., M) 
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47T 

o-totO) = pr Im/pp. (3.6) 

The aim of this section is to derive formulas for the transition and scattering amplitudes for various processes in the 
framework of QED. 

A. Photon emission by an atom 

Consider the process of the photon emission by an atom. According to the standard reduction technique (see, e.g., 
p8| . p6| ), the atomic transition amplitude from state a to state b accompanied by photon emission with momentum 
k/ and polarization €f is 

_i f . ef* exp (ikf ■ y) 
S lf ,b;a = (b\a ou t(kf,ef)\a) = -iZ 3 2 Id y ■{b\j v (y)\a) . (3.7) 

J J2k° f (2w) 3 

Here j v (y) is the electron-positron current operator in the Heisenberg representation, \a) and \b) are the vectors of the 
initial and final states in the Heisenberg representation, Z3 is a renormalization constant, a- b = a v W , e/ = (0, e/), 
kj = (fc?,k/), and k° f = |k/|. Employing the equation (see Appendix A) 

f(y) = exp (iHy°lf(0, y) exp (~iHy°) , (3.8) 

we obtain 



S 7fMa = -iZ z 5 / d A y exp [i(E b + k° f - ^ a )y°]^*(y)(6| > (0,y)| a ) 

= -2mZpd(E b + kj-E a ) fdyA?(y)(b\j„(0,y)\a), (3.9) 



where 



e, exp (ikf ■ x) 
A1(x) = f y{ ' (3.10) 

^2^(2^)3 



is the wave function of the emitted photon. Since \a) and \b) are bound states, equation ( p.9|) as well as the standard 
reduction technique |28Lp6j cannot be used for a direct evaluation of the amplitude. The desired calculation formula 
can be derived within the two-time Green function formalism [g0| |32[ . 

To formulate the method for a general case, we assume that in zeroth approximation the state a belongs to an 
s a -dimensional subspace of unperturbed degenerate states f2 a and the state b belongs to an s^-dimensional subspace 
of unperturbed degenerate states fif,. We denote the projectors onto these subspaces by Pi and P b , respectively. 
We denote the exact states originating from fl a by \n a ) and the exact states originating from fif, by \nb)- We also 
assume that on an intermediate stage of the calculations a non-zero photon mass fj, is introduced. It is considered 
to be larger than the energy splitting of the initial and final states under consideration and much smaller than the 
distance to other levels. 

We introduce 

g^ f {E',E;x' 1 ,...x' N ;yi 1 ,...XL N )8{E' + k -E) 
= — : dx°dx'° / d 4 y exp (iE'x'° - iEx°) exp (ik°y°) 

xAHyKoirvVVi) • • • V-KVW) 

xj v (y)^(x°, x* ) • • • #r°, xi)|0) , (3.11) 



where, as in the previous section, ip(x) is the electron-positron field operator in the Heisenberg representation. Let 

y 7/ in the vicinity of the points E' rs E b 



us investigate the singularities of Q 1; in the vicinity of the points E' rs Ei. ' and E w Ea . Using the transformation 



rules 



:•»:-! 



ip(x°, x) = exp (iHy°)iP{x - y°, x) exp (-iHy°) , 

j(y°, y) = exp (iHy°)j(0, y) exp (-iHy°) , (3.12) 

we obtain 

G lt \E\ E; x' 1: ...xV;xi, ..jf. N )S{Ef + k°-E) 

= / dtdt' / d 4 y exp UE't' - iEt) exp \i(E' + fc° - E)y°] 

2-rri 2tt N\ J _ qo J 

xA v f *(y)(0\TTp(t'y i )---Wy N ) 

xjV(0,y)^(t,x JV )...^(t,xi)|0) 
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xAf(y)(O|T^(f J x' 1 )...-0(t / ,x / w ) 

x>(0,y)^(t,x JV ) • • -?(i,xi)|0) . (3.13) 

Using again the time-shift transformation rules, we obtain 

£ 7/ ( J E / , J E;xi,...x' A r;xi,...x A r) 

= Y~Tn f dtdt ' I d y ex p (^'*' - iEt ) E A r (y) 

x exp (-iE ni t') exp (iE n2 i)0(£')0(-£)<O|2>(O, x'J • • • V(0, x' w )|m) 
x<ni|j*(0,y)|n2)<n 2 |^(0 ) x JV ) • • .^(0,xi)|0) + • • • . (3.14) 

Here we have assumed Eq — 0, as in the previous section. Taking into account the identities 



/ dt exp[i(E'-E ni )t] = 
Jo 
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E' - E ni + iO 



f dt exp [i(-E + E n ,)t] = ' , (3.15) 

J-oo E - hj n , 2 + lU 
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we find 



/ 



4V fdyA V f*(y) — 

N\ ^ J S yy> E' -E ni 



2tt N\ ^ J J J KJ/ E> -Em +ME- E, l2 + iO 

ni,n 2 1 J 



x (0|T^(0, xi) • • • #L x! N )\ni)(ni \j v (0, y) |n 2 ) 

x(n 2 |V(0,Xjv) • ••^(0,x 1 )|0) + • • • . (3.16) 

We are interested in the analytical properties of Q lf as a function of the two complex variables E' and E in the region 

E' « E^ , E w -Eq . These properties can be studied using the double spectral representation of this type of Green 
function which is given in Appendix E (a similar representation was derived in Ref. ]64|). As it follows from the 

spectral representation, the terms which are omitted in equation (3.16) are regular functions of E' or E if E' « E^ 



and E fa Ea ■ The equation ( |3.16| ) and the spectral representation given in Appendix E show that, for a non-zero 
photon mass /j,, the Green function Q 7f (E',E) has isolated poles in the variables E' and E at the points E' = E, lb 
and E = E Ual respectively. Let us now introduce a Green function g~ t ,b-a(E' ', E) by 

g lf , b ; a (E',E) = P^g 7f (E',E)^ ■ ..^?i°) , (3.17) 



where, as in ( p.39| ), the integration over the electron coordinates is implicit. According to equation (3.16) (see also 
Appendix E) the Green function g^, jb . a (E',E) can be written as 

g^ b . a (E',E) = ± £ JT * * <p nb f dyAy{y){n b \UO,y)\n a )ipt a 



n a — 1 rif, — 1 

p(0) 



terms that are regular functions of E' or E ii E' fa E^ 

and J5 « ^ 0) , (3.18) 
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where the vectors tp/. are defined by equation Q2.65 ). Let the contours T a and T(, surround the poles corresponding 
to the initial and final levels, respectively, and keep outside oth er sin gularities of g^ f ^- a ( E / . E) including the cuts 
starting from the lower-lying bound states. Comparing equation (]3.18| ) with equation ( |3.9| ) and taking into account 
the biorthogonality condition ( jj.7l| ) , we obtain the desired formula 
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(3.19) 



where by a we imply one of the initial states and by b one of the final states under consideration. The vectors Vk are 
determined from equations (2.76)- (2.77). 

In the case of a single initial state (a) and a single final state (6), the vectors v a and Vb simply become normalization 
factors. So, for the initial state, 



and, therefore, 
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(3.22) 



we obtain 
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(3.23) 



For practical calculations of the transition amplitude, it is convenient to express the Green function g^, b . a (E',E) 
in terms of the Fourier transform of the 2iV-time Green function, 
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where 
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(3.25) 



The Green function G lf is constructed by perturbation theory after the transition in (3.25) to the interaction repre- 
sentation and using the Wick theorem. The Feynman rules for G 7/ differ from those for G considered in the previous 
section only by the presence of an outgoing photon line that corresponds to the wave function of the emitted photon 
AK* (x) . The Feynman rule for a vertex in which a real photon is emitted remains the same as for a virtual photon 
vertex. The energy variable of the emitted photon (fc°) in the expression corresponding to a real photon vertex is 
considered as a free variable (k ^ k± = |k/|) which, due to the energy conservation, can be expressed in terms of the 
initial (E) and final (E 1 ) atomic energy variables. 
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B. Transition probability in a one-electron atom 

To demonstrate the practical ability of the formalism, in this section we derive formulas for the transition probability 
in a one-electron atom to zeroth and first orders in a. An application of the method for two-electron atoms is considered 



1. Zeroth order approximation 



To zeroth order the transition amplitude is described by the diagram shown in Fig. 31. The formula (3.23) gives 



S\% a = S(E b + kj - E a ) I dE> I dE g^ b . a {E', E) , 



where the superscript indicates the order in a. Here we have taken into account that 



= 1 



and a similar equation exists for the b state. According to the Fcynman rules, we have 

GW((£/x');fc°; (B,x)) = J dy ±S{E' \x\y){-2meY)5{E' + k° - E) 

xA* f>v (y)±S(E,y,x). 
Substituting the expression ( 3.28 ) into the definition of g 1} .b- a {E' , E) (see equation fl3.24| )), we obtain 



g {t) \ (E',E) = - — (b\ea v A* f v \a) — - — 

y lf May ' > 2ttE' -et f ' E-e a 



The equations (|3.26|) , (|3.29| ) yield 



s ylb;a = -2m6(E b + kj- E a ){b\ea v A* ftV \a) 



or, in accordance with the definition (3.2) 



^L = -{b\ea"Al v \a) = (b\ea ■ A}\a) . 
The transition probability to zeroth order is 

dW\% a = 27r|r^ 6;a | 2 <5(^ + k° f E a )dk f 



2tt- 



2fcO(27T) 3 

Integrating over the photon energy, we obtain 



e)*a v exp (-ik/ • *)\a)\ 2 5{E b + k° f - E a )dk f . 



k ( l 



where 



^6 ; a = a^N}-J6a(k/)| 2 dn / , 



J6o(k/) = (b\aexp(-ikf • x) 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



:-!(» 



2. QED corrections of first order in a 



The QED corrections of first order in a are defined by the diagrams shown in Fig. 32. Let us consi der t he derivation 
of the formulas for the self-energy (SE) corrections (the diagrams (a)-(c)) in detail. The formula ( 3. 23] ) gives in the 
order under consideration 

S^, b ; a =S(E b + k° f -E a )U dE' I dEgW b . a (E>,E) 

-\ £ d E> £ dE g^ a (E',E) [± £ dE g^(E) + ± £ dE «&>(*)] } , (3.35) 

where g a <l (E) and g bb , (E) are defined by the first-order self-energy diagram. Here we have omitted a term of first order 

— 1/2 

in a which comes from the factor Z 3 since it has to be combined with the vacuum-polarization (VP) correction. 
Consider first the diagram (a). According to the Feynman rules, we have 



G^((E'yy,k ;(E,K)) = 5(E' + k°-E) / dyd^dz ^-S(E',x',y)^°X(E',y',y 



2tt 



2tt 



x±S(E',y,z)A* fiV (*)(-2merf)±S(E,z,yi) 



2n 



2tt 



where 



£(£', y', y )=e 2 ^-Jdw y°yS(E' - u, y', y)l°D pa {u, y' - y) 



(3.36) 



(3.37) 



is the kernel of the self-energy operator and D pcr {u},y' — y) is the photon propagator for a non-zero photon mass. 
According to the definition of g lf ^ ;a (E', E) (see equation (3.24)), we find 



7/ >b 



«( £ ' £ ) = ^Ew 



(b\Z(E>)\n){n\ea»Al v \a) 



and 



I dE' I dEg ( ^ b ] a (E',E) = -2m[£ 



{E' - e b ){E' - e n ){E - s a ) 

n * b {b\Z{e b )\n){n\e a v Al v \a) 



(3.38) 



where S'(e b ) 



__ rfS(e) 
— de 



£b - £n 

+ (b\Z'(£ b )\b)(b\ea»A}Ja) 
A similar calculation of the diagram (b) gives 

- n * a {b\ea»A* fv \n){n\L(e a )\a) 



(3.39) 



<f dE' I dEg%% a (E',E) = -2iri[£ 



+{b\ea v Al v \a){a\Z'{e a )\a) 



(3.40) 



The s econd ("reducible") terms in equations (3.3E) and ( 3.40| ) have to be combined with the second term in equation 
( |3 .35 ). Taking into account that 



^-£ dEgW(E) = {a\X'(e a )\a), 



2m 



,« 



dEg£>(E) = {b\V( Sb )\b), 



(3.41) 
(3.42) 



we obtain 



dE' 



dE9^, a {E',E)\±-J dEgH\E) + ± 



2iri 



2-ni 



dEg£(E) 



= 2ni 



-(b\ecrA* ftV \a)((b\V(e b )\b) + (a\V(e a )\a)) 



(3.43) 
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For the diagram (c) we find 



9 { 'ijE',E)= ' 



y E ,—j- / rfxdydz Vb(x)— / duj p S(E' -u/,x,z) 



xeA* v (z)YS(E - uj, z, y) 7 <VD p(T ( w , x - y)Va(y) 



1 



£■ 



Substituting ( |3.44|) into ( |3.35|) , we obtain 

I dE' I dE g^f a (E',E) = = -2iri f dzeA* fil/ (z)A»(e b ,e a ,z) 
where 



(3.44) 



(3.45) 



A."(s b ,e a ,z) =e 2 



2tt 



dw / dxdy V' f ,(x)7 P S'(e 6 - w, x, z)7 I/ S'(e - w, z, y) 



x 7 ' T J D pCT (w,x-y)V'a(y). 



(3.46) 



Summing all the first order SE contributions derived above and adding the contribution of the mass-counterterm 
diagrams (Fig. 33), we find 

c(i.se) . x ,„ ,,o tp n \ r ^(b\'Z(e b )-(38m\n)(n\ea v A* f \a) 



^ (&|ea",«»<n|E(e )-/Wm|a) 

+ Z^ Z 

n £a e " 

+ /dzeA}^(z)A y (e b ,e a ,z) 
+i(6|ea"^ iI/ |o)«6|S'( E6 )|6) + (o|E'( £o )|o)) 



(3.47) 



A similar calculation of the vacuum-polarization diagrams (Fig. 32, diagrams (d)-(f)) gives 

^ b (b\Uv P \n}(n\ea»A}Ja) 



s ( ;;;r a ' = -^is(E h +k° f -E a )[j2 



,(1,VP) _ 



£b ~ £r, 



E 



(&|ea v A}»(n|i7vp|a> 



e a -e 



a ^n 



1/2 



where 



is the VP potential and 



+ / dzeAl„(z)Q"(k},z) + (V* - l)(&|eaM^|a 



/-i />00 

rfy, ,/ ^Tr[S( W ,y,y) 7 °] 

Q"(fc°,z) =-e 2 /"dxdy^ 6 (x) 7 Va(x)I? pCT (fc ,x-y) 



(3.48) 



(3.49) 



2vr 



dw Tr[ 7 CT S(w, y, z) 7 "S(w + fc°, z, y)] . 



(3.50) 



Some individual terms in equations ( 3.47 ) and ( 3.48 ) contain ultraviolet divergences. These divergences arise solely 
from the zero- and one-potential terms in the expansion of the electron propagators in powers of the Coulomb potential. 
Using the standard expressions for the divergent parts of the zero- and one-potential SE terms (see, e.g., |)28|) and 
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the Ward identity (Z\ = Z2) one easily finds that the ultraviolet divergences cancel each other in equation (3.47). As 



to equation ( |3.48| ), the divergent parts incorporate into the charge renormalization factor (e = Z 3 ' eo) 



An alternative approach to the renormalization problem consists in using the renormalized field operators 



if'R 



- y-1/2. 



- 5T-1/2 



tp, An — Z 3 A, the renormalized electron charge e = eo + Se = Z x Z<xZ 3 eo, and, respectively. 



1/2, 



the renormalized Green functions from the very beginning. In that approach, additional countcrterms arise for the 
Feynman rules. 



The vertex and reducible contributions to the SE corrections (third and fourth terms in equation (3.47)) contain 
infrared divergences which cancel each other when being considered together. 

In addition to the QED corrections derived in this subsection, we must take into a ccount the contribution originating 
from changing the photon energy in the zeroth-order transition probability (3.33) due to the QED correction to the 
energies of the bound states a and b. It follows that the total QED correction to the transition probability of first 
order in a is given by 



dW, 



(i) 



7/, 



2^) 2 2Re{r^V« b;Q }^ / + 



dW, 



(0) 



7/, 



k°=E a 



dW r 



(0) 



ff< 



k f -e 



(3.51) 



(i) 

jf ,b;a 



(1,SE) 

~/f,b;a 



(l.VP) 
7/,6;a 



is the QED correction given by equations (|3.47| ) and ( |3.4<j ) in accordance with the 



Here r, 

definition (3.2). E a , Ef, and e a , Ef, are the energies of the bound states a and b with and without the QED correction, 

respectively. 



C. Radiative recombination of an electron with an atom 



In calculations of processes which contain a free electron in the initial or final or in both states we consider that 
the interaction with the Coulomb field of the nucleus Vc(x) is included in the source j(x) which leads to a scattering 



(i $ -m)ip(x) = j(x) . 



(3.52) 



It means that, after the transition to the "in" operators, the unperturbed Hamiltonian does not contain the interaction 
with the Coulomb potential. As a result, the Feynman rules contain free-electron propagators, instead of the bound- 
electron ones, and the vertices corresponding to the interaction of electrons with Vc(x) appear. Since we consider the 
case of a strong Coulomb field, we will sum up infinite sequences of Feynman diagrams describing the interaction of 
electrons with the Coulomb potential. As a result of this summation, the free-electron propagators are replaced by 
the bound-electron propagators, 



[p° - H{1 ^O)]- 1 = [p° - H (l - tO)}- 1 

+\p°-Ho(l-iO)]- 1 Vc\ P ° 



H o (l-i0)] 



(3.53) 



where H = Hq + Vc and Hq — —ia ■ V + /3m, and the free-electron wave functions are replaced by the wave functions 
in the Coulomb field. For instance, the wave function of an incident electron with momentum p^ and polarization \Xi 
is 



V> P , 



M+) 



U, 







¥1 



u v 



- {p° - H (l - iO^VcU^ 

#0(1 - *o)]-V c [p° - #0(1 - zo)]- 1 ^, + ■ 

-[p?-#(l-*0)]- 1 F c CV J , 



where pi = (p°,Pi 



vw 



u. 



Pi m 



M(p,,/i t )exp(ipi -x) 



(3.54) 



(3.55) 



is the free wave function of the incident electron, and u(j>i, m) is normalized by the condition uu = 1. 

Consider the process of the radiative recombination of an electron with momentum p^ and polarization fii with an 
(N — l)-clectron atom X^ Z ~ N+1 ^ + in a state a. As a result of the process, the N-electron atom X( Z ~ N ^ + in a state 
b and a photon with momentum k/ and polarization e/ = (0, e/) arise, 



"(Pi,*) + X( z - N +V+(a) - J(k f ,e f ) + *<*-">+(&) ; 



(3.56) 
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where kf = (fc9,k/) and k° f — |k/|. In this section we will assume that we consider a non-resonant process. It means 
that the total initial energy (p° + E a ) is not close to any excited-state energy of the TV-electron atom. As to the 
resonance recombination, a detailed description of this process is given in pT]] (see also the next sections of the present 
paper). 

According to the standard reduction technique Pq,p6| , the amplitude of the process is 

S lf ,b- ei ,a = (b\a out (kf,ef)bl n (pi,Hi)\a) 

ef exp (ik f ■ y) 



= (-iZ 3 5 )HZ 2 5 ) / d 4 yd 4 z € 



2k° f (2n) 



x (b\Tj v {y)ip{z)\a){-i Z -to) . = , (3.57) 



where \a), \b) are the vectors of the initial and final states in the Heisenberg representation. Taking into account that 

V>(z)H *0z - m ) = e${z)A(z) = rj(z) (3.58) 



we obtain 



5 7/>6;ei , a = 2n5{E b + k° f - E a -p° i )(-iZ 3 h )(-iZ 2 h ) J dydz Af(y) 
x{ / dr (b\j u (T,y)r](0,z)\a) exp (ik° f T) 



dr (%(r, z) j„ (0, y) | a) exp (-ip°r) 



ii 



t(6|[^(0,a), ju(P,y)]\a)}U PilH (m) , (3.59) 



where A v Ay) is the wave function of the emitted photon defined by equation ( 3.1 0| ) and U PifJii (z) is the free wave 



function of the incident electron defined by equation ([3.55). To construct the perturbation theory for the amplitude 



of the process, we introduce the Fourier transform of the corresponding two-time Green function 

G 1} , et {E\ E,p°; x' 1; „.xk;xi, ...x N -i)8(E' + fc° - E - p°) 

= (±V I l f°° dx°dx'° exp(iE'x'°-iEx 

\2ir) y/Niy/iN-iy.J-co 

j4„.j4_ /".'I.0..0 ■ 0_0\ ai/* 



x / dTyefz exp(ik v y u -ip v z v )A'j*(y) 

x(0|T^(x'°,xf) • • •^(x'°,x , w )^(y)^(z)^(x°,x Ar _ 1 ) • • ■^(x°,^ 1 )\0) 
x(-if)z-m)U pm {z). (3.60) 

As in the derivation of the formulas for the transition amplitudes, we will assume that in zeroth approximation the 
initial and final states are degenerate in energy with some other states and we will use the same notations for these 
states as in section III(A). As usual, we also assume that a non-zero photon mass [i is introduced. The spectral 
representation of G^ f ; ei (E' , E,p°) can be derived in the same way as for the Green function describing the transition 
amplitude (see section III(A) and Appendix E). It gives 

g if . ei (E',E,p°;-x' 1 , ...x'^jxi, ...xjv-i) 



1 1 



== y y — - 



2vr VM y/(N-l)\ ^ ^ E> - E nb E - E na 
x fdydzAf{y)(0\TTp{Q,x!?)---i>{Q,x! N )\n b ) 



x{ / dr (n b \j u (T,y)r](0,z)\n a ) exp{i(E - E nb + p u )r) 
Jo 
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+ dr (n b \i](r,z)j u {0,y)\n a ) exp{i(E' - E m - p°)r) 



-i(n b \[^(0,z), j v (0,y)}\n a )}(n a \^(O^N-i) ■ ■ ^(O^MUp^i 
+ terms that are regular functions of E' or E if E' sw E^ 
and E w ^ 0) . 

We introduce the Green function g(£", E,p°) by 

ff7/ , 6;es , a (£',£:,p ) = PJ; 0) g^ i e i (E' 1 E,p )^ ■ ■•7^-i^ (0) • 
^From equation (3.6l|) we have 



(3.61) 



(3.62) 



g^ b . ei , a (E',E,p°) = l- £ JT E ^ b E jty^Ayb) 



n a — 1 rib — 1 

/•OO 

x{ / dr (n 6 |j 1 ,(T,y)r/(0,z)|7i a )exp(i(£ , -£ , nb +p°)r) 
dr (n 6 |^(r,z)>(0,y)|n a )exp(«(i5'- i5 ria - p°)r) 



-z(n 6 |[^ t (0,z),^(0,y)]|n a )}[/ PiAl! (z) 



_E — EL 



terms that are regular functions of E' or E ii E' m E^ ' 
and £ « £i 0) , 



(3.63) 



where ip Ua and ip nb are the wave functions of the initial and final states as defined above. Let the contours T a and 
F b surround the poles corresponding to the initial and final levels, respectively, and keep outside other singularities of 
g-y f ,b;ei,a(E', E,p°). This can be perfor med if the photon mas s ji i s chose n as i ndicated in section III(A). Taking into 
account the biorthogonality condition ( 2.71 ) and comparing ( 3.63 ) with ( |3.59 ) we obtain the desired formula [B7J 



Sj f , b ;e i ,a = (Z2Z 3 )-i8(E b + k f -E a -p^f dE' A dEvlg^ t , b . ei ,a{E',E^)v a , 



(3.64) 



where by a we imply one of the initial states and by b one of the final states under consideration. In the case of a 
single initial state (a) and a single final state (b), it yields 

S JJ , b -e l ,a^(Z 2 Z 3 )-h(E b + k f -E a -p t )(f dE 1 I dEg lsM ^ a (E\E lV \) 



-^ d> dEg bb (E) -^ <£ dEg aa {E) 



2m 



2m 



-1/2 



(3.65) 



For practical calculations by perturbation theory, it is convenient to express the Green function <7 7/ ,6 :ei ,a in terms 
of the Fourier transform of the (2N — l)-time Green function, 

g^, b . ei ,a{E', E,p°)S(E' + k°-E-p°) 

i r°° 

===== / dp° v -- dp%_ x dp>? ... d p%5{E-p\---- p%_ x )5{E> -pf---- P %) 



^/n\(n-i)\ 

b W/; 



yP (0 V fV° r,'°-k° rfi-rfi t>° W ...^ p(°) 



(3.66) 



where 



G 1 , e ,((pr,x' 1 ),..,(|,S,x^);iV ;(P^i),.- I (l'li ! XN-i)) 






— ) (271)2^+! J_ OQ dX °'" dx °N-l dx 'l ■■■dx'l 

x exp (ipfx? + ■■■ + ip%x% - ip ^ ip%-i^N-i) 
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x I d^yd^z exp(ik°y° - ip°z°) Af{y) 

x(-i0 z -m)U pulz {z). (3.67) 



The Green function G 7 , ei is constructed using the Wick theorem after the transition in (3.67) to the interaction 
representation. Since we have not included the interaction with Vc in the unperturbed Hamiltonian, the Feynman 
rules contain the free electron propagators and the vertices corresponding to the interaction with Vq- Summing over 
all the insertions of the vertices with Vc in the electron lines we replace the free electron propagators and wave 
f unct ions by the electron propagators and wave functions in the Coulomb field, according to equations ( 3.53| ) and 



(3.54). The free wave function of the incident electron U PifJii (x) is replaced by the wave function ip Pil j, i (+){^) that can 



be determined from the equation 

vw*(+) = Up** + (p? - w - *o)rv c ^M+) ■ ( 3 - 68 ) 

The wave function tp ViUi (+)(' x ) is a solution of the Dirac equation with the Coulomb potential Vc(x) which satisfies 
the boundary condition 81 



Vw( +) (x) ~ LWx)+G+(n) CXp( ff |x|) , M-oo, (3.69) 

where G + (n) is a bispinor depending on n = x/|x|. The apparent expressions for ^pu(+)(x) are given, e.g., in 
p2| . Thus, we again revert to the Furry picture. The Feynman rules for G 7/;ei differ from those for G 7/ only by a 
replacement of one of the incoming electron propagators ^-S(u>,x, y) by the wave function ^ p . M .( + )(x). 

In calculations by perturbation theory a problem appears which is caused by the fact that in zeroth approximation 
the energy of the (TV— l)-electron atom may coincide with the difference o f the iV-electron energy and the one-electron 



energy. As a result, some of the terms which are omitted in equation (3.61) are singular functions of E' and E if 



E' rs E\ and E ra E a . A special analysis of the complete spectral representation of the Green function Q lf , ei shows 



that to eliminate these terms in calculations by formula (|3.64) the following prescription should be used: integrate 



first over E and keep the point E = Ea + (£" — E^ ') outside the contour T a . 

In Refs. p3,p6[ this method was used to derive formulas for the QED corrections to the radiative recombination of 
an electron with a bare nucleus and for the interelectronic-interaction corrections to the radiative recombination of an 
electron with a high-Z heliumlike atom. These formulas are presented in sections IV(D1), IV(D2). Here we consider 
another application of the method. 

D. Radiative recombination of an electron with a high-Z hydrogenlike atom 

As a practical application of the method, in this section we derive formulas for the radiative recombination of an 
electron with a high-Z hydrogenlike ion to zeroth and first orders in l/Z. We will assume that the final state of the 
heliumlike ion is described by a one-determinant wave function 

u 6 (xi,x 2 ) = — ^(-l) P V>P6 1 (xi)^Pb 2 (x 2 ) (3.70) 

and the one-electron state |6i) coincides with the initial state \a) of the hydrogenlike ion. 

1. Zeroth order approximation 
To zeroth order, the radiative recombination amplitude is described by the diagram shown in Fig. 34. The formula 



(3.65) gives 

s!$ Mei ,a ^ 5 (^ +k° f -E a - p°) I dE' I dE g^ >b , ei>a (E', E, P °) , (3.71) 



where the superscript indicates the order in a. According to the definition ( 3.66| ) and the Feynman rules, we have 
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1 J>lf r dp[°dp' 2 ° 8{E' - p? - p' 2 °) ' 



2tt 



x {Pb 2 \ea u A v f * \p)5(p' 2 v + k u - pY 



1 



pf -e a + iQ 



p' 2 - £pb 2 + iO 
Sip? - E){Ph\a) 



where \p) = \pi, Hi) denotes the state vector of the incident electron. We obtain 



C b ^J E '> E ^ 



1 



■(b 2 \ea u Af\p)- 



1 



2ttE' -E-e b2 N ' J > r 'E-e a 

Following to the rule for the integration over E and E 1 formulated at the end of section III(C), we find 



<f dE' I dEg\ f \ b . et jE\E 7 p° l ) = -2m(b 2 \ea u AY\p) 



? (0) 

7/,6;ei,a 



It yields 



or, according to the definition (|3.2|), 



The corresponding cross section is 



-2iri5(E b + k° f ~ E a - p° t ){b 2 \ea v Ay\p) 



(0) 

■jf,b;e 



■{bi\ea v A v f * \p) . 



daW (2tt) 4 (0) 2 

where Vi is the velocity of the incident electron in the frame of the nucleus. 



(3.72) 



(3.73) 



(3.74) 



(3.75) 



(3.76) 



(3.77) 



2. Interelectronic-interaction corrections of first order in 1/Z 

The interelectronic-interaction correction s of fi rst order in l/Z arc defined by the diagrams shown in Fig. 35. In 
the order under consideration, the formula (3.65) gives 



S^ Mi , a = S(E b + k° f -E a -p9) j> dE' f dEgW b . eita (E',E,p°) 



i 



i 



dE' f dEg ( °l b . !ei jE\E,p°) — f dEgg\E) 



(3.78) 



where g bb (E) is defined by the first-order interelectronic-interaction diagram (see Fig. 20). According to the definition 
(|3.66) and the Feynman rules, we have 



9^, b ., ei .a^',E,p^5{E' + fc° - E p°) 

= (£) E(-!) p J *$*£ S ( E ' pi - p'") 



i 



i 



i 



x 



dq • duj 



Pi - e Pbl + iO p' 2 ° - s Pb2 + iO E - e a + iO 
x J2 [(PbiPb 2 \I{u;)\an)5{p 2 + u; - q°)5(pf - u - E) 

n 

l 



{n\ea v A)*\p)5{q v + k»-p») 



q°-e n {l-i0) 
+ (Pb 1 Pb 2 \I(u;)\np)5(p' 2 +00- p°)8(ft? - w - q°) 
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1 



g°-e„(l-tO) 

-{Pb 2 \ea v Ay\n) — 



{n\ea v Ay\a)5{q a + k° - E) 
1 



-Ap 



10 ' k° - q° 



We obtain 



<T -s„(l-*0) 

x (Pb in \I(u;)\ap)S(q +LU- p1)5(p'? - u - E) 

+(Ph\ea v A?\n) q0 _^ {1 _. 0) 6(p? + k° q°) 

x(nPb 2 \I(Lj)\ap)5(p'° + u - p°)S(q° - uj - E) 



■ 9 1 1 poo 

«™w<*.*..*> = (£) Fr^r^:? ( - 1) '?/-* f 



(3.79) 



Pi u - e Pbl +t0 E< - p'l> - e Pb2 + »0 

(PhPhim - e)\c • ' 



E'-E-s n (l-i 
+ (Pb 2 1 ea u A v f * I n) R «i ; r 

x(P6m|/(p / 1 °-£)|ap) 

/ i N2 



— (nlea^A^lp) 



■ - 



V27r/ £'_£ 



^^E(-^E/> 



f l 



V 2 - £pb 2 +i0 E>- p'$ - e Pbl + tO 



(P6 1 P6 2 |/(rf-p' 2 u )|np)- 



1 



-(PftileavA'fln) 



£' - pi - e n (l - iO) 
1 



(n\ea v Ay\a) 



f ry E + p°-p> 2 °-e n {l-iQ) 
x(nPfe 2 |/(pO-p 2 °)|ap)" . 



(3.80) 



The first term in equation ( |3.80|) is conveniently divided into an irreducible (e n ^ E b ' — e a = £& 2 ) and a reducible 
( e n = Eb 2 ) part. Since we consider the case of a single final state described by the one-determinant wave fu nction 
( 3.70) wit h |&i) = \a), the condition e n — e b2 means \n) = \b 2 ) (otherwise (PbiPb 2 \I(uj)\an) = 0). Substituting ( 3.80J ) 
into ( 3.78 1 ), int egrating over E and E' according to the rule formulated at the end of section III(C), and using the 
identity (p. 127), we find for the irreducible contribution 



S%^ = 6(E b + k° f -E a -p°)f dE> f dEg^Z(E',E,p°) 



(l,irred) , „, 



= -2mS(E b + fc° - E a -p°) J2(- l Y 



n=£b 2 



[j2{Pb 1 Pb 2 \I{£p bl -e a )\ 

n 

-J2(Pb2\ea u A1 



£bo £-r, 



(n\ea v Ay\p) 



■ — n (Pbm\I(ep bl - £ a )\ap) 

£a +Pi —SPb! —£n 



J2{PbiPb 2 \I(p t - e Pb2 )\np) l (n\ea u Af\a) 



& b Pi tn 



■Y^{Pb x \e aiJ Ay\n) - l 

*-* J £a+Pi ~ £Pb 2 - £r. 



(nP6 2 |/(p°- ep62 )M 



(3.81) 



For the reducible part we have 



44 



f dE> I dEg^Xl a {E',E^) 
dE> £(-1 



T, 

i \2 



{E'-E^f 



dp? 



1 



Pi u - s Pbl + id 



E'-p?-e Pb2 +iO 



(Pb 1 Pb 2 \I(p? - e a )\ab 2 ){b 2 \ea v A?\p) 



■£(-i) p f 

D J —I 



du) 



(w - epb! - i 



(P6 1 P6 2 |J(a; - e a )\b x b 2 ){b 2 \ea v Af\p) . 

1 — i\)) z j 



(3.82) 



The reducible contribution must be considered together with the second term in formula ( p. 78 ). Taking into account 
that 



-L/ d Eg${E)=-i- 



2iri 



2- 



2 / dp? 

J — OO 



1 



(p'° - e bl - zO) 2 
1 



(6i6 2 |/(K° - e6j|6i6a> 



(p? - e b2 - zO) 2 



(b 2 bi\I(p? - e bl )\bib 2 ) 



we find 



:i(W«uy' W r.r i .(h«»W»)l».W 



(3.83) 



do; 



(w - z'O) 2 

(6 a 6 1 |/(w)|6i6a)( 7 * ^ + T , , ..., 

V (uj — A b — z0) z (w + Ah — i0) z 



1 



(3.84) 



where Af, = s b2 — e bl - Summing (3.82) and (3.84Q, we obtain for the total reducible contribution 



-,(l,red) _ 
7/ : b;ei ,a 



1 7 f°° 

-2iri6(E b + k f-E a -p1)~{b2\ea v A?]p) — J dto (& 2 &i|/(u#i& 2 ) 



(w + A 6 + iO) 2 (w + A 6 - iO) 2 



Using the identity 



(w + iO) 2 {uj - iQ) 2 



2ir d 
i duj 



8{lu) 



and integrating by parts, we obtain 



S< iXla = 2mS ( E b + k° f - E a -p^-(b 2 \ea,Ay\p)(b 2 b 1 \I'(A b )\b 1 b 2 ) 



(3.85) 



(3.86) 



(3.87) 



In addition to the interelectronic-interaction correction derived in this subsection, we must t ake i nto account 
the contribution originating from changing the photon energy in the zeroth-order cross section ( p. 77 ) due to the 
interelectronic-interaction correction to the energy of the bound state b. It follows that the total interelectronic- 
interaction correction to the cross section of first order in 1/Z is given by 



da™ (2tt) 



dflf 

(1) _ a.iri'cd) (l,red) 



•II k 2 2 Re(T (0) * r (1) X 



rdo-W 



<m 



da^ 



fc°=p° +£„-£,, dfi 



f 



k o p o +Sa _ E (o) 



(3.88) 



Here t:, ' b . e a = r^ 'g * „ + T2. '™ e a is the interelectronic-interaction correction given by equations ( |3.8l|) and ( |3.87j ) 



in accordance with the definition (I3J). E b and St are the energies of the bound state b with and without the 



interelectronic-interaction correction, respectively. 
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E. Photon scattering by an atom 

Let us consider the scattering of a photon with momentum k^ and polarization e, = (0, €{) by an atom which is 
initially in a state a. As a result of the scattering, a photon with momentum kj and polarization ey = (0, €/) is 
emitted and the atom is left in a state b. In this section we consider the non-resonant photon scattering. This means 
that the total initial energy of the system (E a + kf) is not close to any excited-state energy of the atom. The resonant 
photon scattering will be considered in detail in the next section. 

According to the standard reduction technique (see, e.g., p8fl ), the amplitude of the process is 

S-y f> b- ni , a = {b\aont(kf-,ef)ai n (ki,ei)\a) 

_1 f 4 4 e f ex P (*^/ ' V) 

= Disconnected term — Z 3 Id yd z — — . — 

J J2k f {2TT) 3 

x(b\Tj„(y)j p (z)\a) e *™l { -^- z) , ^ 

V 2k i( 27r ) 

where \a) and \b) are the vectors of the initial and final states in the Heisenberg representation; kf — (fc?,ky) and 



ki = (kf, hi). The first term on the right-hand side of equation (3.89) corresponds to a non-scattering process, i.e., the 
photon does not interact with the atom. We are interested in tr. 



by the atom. Let us designate this term by S^ cat 



le second term which describes the photon scattering 



To derive the calculation formula for Sff^l a , in the scattering amplitude we isolate a term which describes the 
scattering of the photon by the Coulomb potential Vc- With this in mind, we write 

(b\Tj(y)j(z)\a) = (b\[Tj(y)j(z) - (0\Tj(y)j(z)\0)]\a) +6 ab (0\Tj(y)j(z)\0) . (3.90) 

The second term on the right-hand side of this equation corresponds to the photon scattering by the Coulomb field 
and the first term describes the photon scattering by the electrons of the atom. To the first term, only diagrams 
contribute where the incident photon is connected to at least one external electron line. We denote this term by 
S c ™ b ^ a - We have 

-i f ^ a e/exp(%-y) 



S5$. 7li0 = -Z3- 1 / d*yd 4 z J V J ' (b\[Tj v (y)j p (z) 

/2fc°(2^) 3 



-(0\Tj„(y)j p (z)\0)}\a) € * eX llj t=^ 
\ 1 Jv\y>j P \ n /ji / ^ 2jfc o( 2 ^)3 

/OO P 

dt exp(ik° f t) / dydz Af(y) 
-OO J 

x(b\[Tj u (t,y)j p (0,z) - (0\Tj„(t,y) Jp (0,z)\0)}\a)AP(z) . (3.91) 

Here the second term is zero if a 7^ b. 

To construct the perturbation theory for the amplitude of the process, we introduce the Fourier transform of the 
corresponding two-time Green function 

S 7 ™ ,(£', E, k'°; x' 1; ...x'^Xi, ...x N )6(E' + k'° - E - k°) 

■ o 1 POO p 

= (_L j _£. / dx dx >0 / d iy d i Z exp (itfgjQ + ^'0^0 _ ^^ _ . fc z0) 

xAf(y)A$(z)(0\TTP(x'\x' l )--^(x'°,x' N )\j v (y)j p (z) 

-(0\Tj,(y)j p (z)\0mx°^ N ) ■ ■ .^(x°,Xi)|0) . (3.92) 

Introducing t 1 = x — z , t = x° — z° , and r = y° — z° and integrating over z°, we obtain 

g^ yi (E',E, k'°;x[, ...x^; Xl , ...k n )6(E' + k'° - E - k°) 

1 1 f°° 

= -8{E' + k'°-E-k°) / dtdt'dr exp (iE't' + ik'°T - iEr) 

2tt N\ J_ oc 

xjdydz Af (y)^(z)(0|T^(i',xi) • ■ • V(i',x^)[^(r, y ) Jp (0, z) 

-(0|Tj v (t, y)j P (0, z)|0)]^(t, XA r) • • • ^(i, Xl )|0) . (3.93) 

4G 



Integrating over t and t' we find 



g™ li (E',E,k'°;x' l ,...x' N ;x u ...yL N ) 



1 1 

2nN\^ E'-E n + iOE- E m + iO 
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/oo 
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-oo 
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■<0|Ti v (T,y)i P (0 



PV"j- 



m)(m|V(0,x^)---V(0,x 1 )|0) + •••}. 



(3.94) 



It can be shown that the terms which are omitted in equation (3.94) are regular functions of E' or E if E' « E'b 
and _E w _E a . As in the previous sections of this paper, we assume that in the zeroth approximation the initial and 
final states of the atom are degenerate in energy with some other states and use the same notations for these states 
as above. As usual, we also assume that a non-zero photon mass (i is introduced. We define the Green function 



57/ , 6;7 „ Q (£', E, k>°) = P^g™ i (E', E, k'°H ■ ■ ■ 7 ^i 0) 



(3.95) 



From equation (B.94) we have 

g^ b .^ a (E',E,k'°) = — J2 E ^ h E j d y dz AT(y)A p l (z)j dr exp (ik'°r) 



x{9(T)e X p[iT(E' - E nb )](n b \\j v (T,y)j p (0,z) 

-(0\Tj u (T,y) Jp (0,z)\0)}\n a ) + 9(-T)exp[ir(E na - E)} 



x(n 6 |[j p (0,z) > (r,y) - (0\Tj v ( T ,y)j p (0,z)\0)}\n a )} 



<P 



E — E n 



terms that are regular functions of E' or E if E' sw E- { 
and E « E { a ] . 



(0) 



(3.96) 



Taking into account the biorthogonality condition ( |2.7l| ) and comparing (3.91) with (3.96), we obtain the desired 
formula |3Q1 



fi^. 7J> „ - Z^S(E b + k°, -E a - *°) <f dE' I dE vlg^jE 1 , E, k° f )v a 



(3.97) 



where we imply by a one of the initial states and by b one of the final states under consideration. In the case of a 
single initial state (a) and a single final state (6) it yields 

S^ biJi , a = Z^5{E b + k° f -E a - fc°) I dE' I dE 9 ;° n bnt jE\ E, k° f ) 



1 



dEg bb (E) 



-1/2 



1 

2~7ri 



dEg aa {E) 



-1/2 



(3.98) 



The disconnected term describing the scattering of the photon by the Coulomb field is calculated by the formula 

(3.99) 



odiscon 

^7_f ,6;7i,a 



-Z 3 1 $ab I d A yd 4 z 



e ^ p{lkr %\TMy)M^ e ^- lh - z) 



^kj(2nr 



vm^r 



For practical calculations by perturbation theory it is convenient to express the Green function g 
of the Fourier transform of the 2A-time Green function, 



7 ° n &; 7 „a in termS 
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C^,.^, E, k'°)S(E' + k'° E k°) 
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v! , dpi--- dp^dp'? ...dp'ZSiE-p*---- p° N )5(E> -p[°---~ p#) 
xP^G^M , -;Pt k'°, fc°;rf, ...,^) 7l ° • • • 7 ^Pi°) , (3.100) 



where 
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x exp (ipfxf + ■■■ + ip%x% - ip\x\ ip%x° N ) 

x f d 4 yd A z ex-p{ik' a y°-ik a z a ) A'/iy) 

x(0|TV-K) • ■■i>(x' N )\j v (y)j p (z) - (0\Tj v {y)j p (z)\0)] 

x^(x N )---^{x 1 )\0)AP{z). (3.101) 



The Green function G 7 on 7 . is constructed using the Wick theorem after the transition in ( 3.101 ) to the interaction 
representation. The Feynman rules for GS° n differ from those for G 7 , only by the presence of the incoming photon 
line which corresponds to the incident photon wave function A?(x). 

F. Resonance scattering: Spectral line shape 

The formulas for the scattering amplitudes derived above allow one to perform calculations by perturbation theory 
in the case where the total initial energy of the system is not close to the energy of an intermediate quasistationary 
state. This is the so-called non-resonant scattering. In the case of resonance scattering, when the initial energy of 
the system is close to an intermediate-state energy, the direct calculation by perturbation theory according to the 
formulas derived above leads to some singularities in the scattering amplitude. It means that those formulas cannot 
be directly applied to the resonance processes. In this section we formulate a method which allows one to calculate 
the resonance-scattering amplitudes. This method was worked out in Ref. p3| . Another approach, which is limited 
to the one-electron atom case, was previously developed in Ref. |p3| . An attempt to describe a decay process within 
QED was undertaken in Ref. J84J. 

The photon scattering by an atom that is initially in its ground state a is now considered for the case of resonance 
E a + k® ~ Ed(d = 1, ..., s), where E a is the ground state energy of the atom, kf is the incident photon energy, and 
Ed(d = 1, ..., s) are the energies of intermediate atomic states which in zeroth order approximation are equal to the 

unperturbed energy E\ ' of a degenerate level. We consider that, as a result of the scattering, a photon of energy 
kf = k® is emitted and the atom returns to its ground state a. The calculation of the photon scattering amplitude by 
the formula derived above leads to a singularity which is caused by the fact that in any finite order of perturbation 
theory one of the energy denominators of an intermediate Green function is close to zero. Therefore, in the calculation 
of the Green function g™ n a (E', E, k'°) we have to go beyond the finite-order approximation. With this in mind, 

let us represent ff^™ 0;7lia (-E', E, k'°) as 

ff ™ n a;7i! a (E', E, k'°) = ^g a (E')Rtf(E', k'°, E + k°)^g d (E + k°)R^(E + k a ,k a ,E)^-g a (E) 

+Ag^ a .^ a (E',E,k' a ), (3.102) 



where g a are gd are the Green functions defined by equation (p. 63) with the projectors P a and P^ , respectively, 



k° = k'° + E' - E, and Ag™* a . lua (E' , E, k'°) is a part of the Green function ff 7 ° n a:7iiQ (-B', E, k'°) which is regular at 
E + k a ~ Ed(d = 1, ...,$). The operators R\ s and i? 7i are constructed by pe rturb ation theory from equation (3T02) 



which must be considered as their definition. Taking into account equation ( p. 64 ) und using the formula ( 3.97 ) , we 
obtain 

S™ n Q;7l , a = V*(*° " k?)cplB!-\E a , k% E a + k°)^-g d (E a + kf)R^(E a + fc°, fe°, E a )<p a 

+Z- l 5{k°-kf)j dE' j dEvlAg™ a . lua (E',E,k° f )v a . (3.103) 
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Consider now how the intermediate Green function gd(E a + k°) can be calculated. Let us introduce a quasipotential 
V d (E)by 



,((» 



(o). 



g d (E)=g d »>(E) + g d »>(E)V d (E)g d (E), 



(3.104) 



where g d = P d /(E— E d ). The quasipotential V d (E) is constructed by perturbation theory according to equation 



( p.l04|) which must be considered as its definition. This equation yields 

V d (E) = [g d 0) (E)]-i - ME)]' 1 = [^(E)}- 1 - [gf{E)+g d 1 \E) 

= [9 d ° ) (E)rg d 1 \E)[ 9d a \E)r + - 



(3.105) 



If the quasipotential V d (E) is constructed from equation ( 3.105| ) to a finite order of perturbation theory, the Green 
function g d is determined by 



g d (E) = [E- Ef - V d {E) 



(3.106) 



The Green function g d (E) has poles on the second sheet of the Riemann surface, slightly below the right-hand real 
semiaxis (see Fig. 5), and has no singularities for real E when E ~ E d . It means, in particular, that if we take the 

quasipote ntial a t least to the lowest order approximation ( V(E) « V(E d ) ), the Green function g d (E) calculated by 
equation (3.106) has no singularities at E ~ E d {d — 1, ..., s), and, therefore, the calculation of the resonance-scattering 
amplitude by equation ( |3.103j ) will be correct. The calculation of g d (E) by equation (3. 10£) effectively corresponds 
to summing an infinite subsequence of Feynman diagrams. 

For the calculation of g d (E) to the lowest order approximation it is convenient to introduce an operator H. by 



H = E 



(0) 



V d (E d ' 



(oh 



(3.107) 



The operator 7i is not Hermitian and has complex eigenvalues. We assume that Ti is a simple matrix, i.e., its 
eigenvectors form a complete basis in the space of the unperturbed d-states. We denote its eigenvalues by £ d = 
E d — iT d /2, the right eigenvectors by \d R ), and the left eigenvectors by |dz,). It means 



H\d R ) = £ d \d R ) , {d L \H = {d L \£ d . 

It is convenient to normalize the vectors \d R ), \di) by the condition 



(d L \d' R 



Tdd' 



(3.108) 



(3.109) 



They satisfy the completeness condition 



d=l 



\d R )(d L 



For g d (E) we obtain 



g d (E)*(E-H)-i = J2^ 



d =i 



\ dn){d L \ 
E-S d 



(3.110) 



(3.111) 



In fact, due to T-invariance, Tiik — Tiki- For this reason the components of the vector (d^l can be chosen to be equal 
to the corresponding components of the vector \d R ). In other words, the components of the vector \d R ) are equal to 
the complex conjugated components of the vector \di). If the d states have different quantum numbers, such as the 
total angular momentum or the parity, one finds \d R ) — \di) = \d). 

Substituting the lowest order approximation for g d (E) given by equation (3.111) into ( 3.103| ), we obtain in the 
resonance approximation 



crcon 



I 

2^ 



6(k° f - fc 4 °) 



o\ Y^ 



fci Ea 



k? -E d + iT d /2 



(3.112) 



In the resonance approximation, it is sufficient to eval uate th e operators i? 7/ , i? 7i to the lowest order of perturbation 



theory. They are determined directly from equation (3.102). 

To demonstrate how the method works we consider below the resonance photon scattering on a one-electron atom. 
A more general case of a few-electron atom is considered in Q . 
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G. Resonance photon scattering by a one-electron atom 

To study the resonance photon scattering by a one-electron atom, in the lowest order we must consider the diagram 
shown in Fig. 36. The contribution of this diagram to the Green function G 7/i7i is 

dydz / dp —S(E',^,y) — e 1 J(E' + k /0 -p°) 



' — oo 
i „, n , 27T 



xAY(y) — S(p ,y,z) — e lp S(p° - k° - E) 

1 Z7T % 



x^(z)^-S(£,z,x). (3.113) 

For ff7°?a ;7l ,a w e obtain 

± WA Jl a p±MK . (3.114) 

2tt ^ E + k° - e n (l - iO) i p l 2tt E - e a K ' 

We consider that E a + k® ~ Ed and, therefore, represent g™ n a . a (E r , E, k /0 ) as the sum of two terms 
Cn,.( E '^ fc '°) = ^\E') 2 -1 ea^Afi-g^iE + k ) 

x if .^o'^L^ ^^^) • (3-115) 



7). 



£; + fc -£„(i-iO) 



Comparing this equation with equation ( 3.102 ), we derive 



9tt 9tt 

47 } (y) = -jea v Af (y) , i? 7 +>(z) = — ea p Af(z) . (3.116) 

Let us derive now the quasipotential Vd(E). To the lowest order of perturbation theory it is defined by the SE and 
VP diagrams (see Figs. 12,13). As was derived above (see section 11(E)), the contribution of these diagrams to ga{E) 
is 

ti\E) = ^ 0) (£ SE (25) + U^gf. (3.117) 

Therefore, 

V} 1] (E) = [$\E)]^g { l\E)[g { P{E)]^ = P^ (X SE (E) + U VP )P^ . (3.118) 

The operator P d Y>s~e(E)P^ contains a non-Hermitian part which is responsible for the imaginary part of the energy. 
The operator 

H = e d + V d [1 \e a + k°) (3.119) 

acts in the s-dimensional space of the unperturbed states. In reality, due to the fact that the operators Sse and f/yp 
do not mix states with different quantum numbers and among the degenerate one-electron states there are no states 
with the same quantum numbers, in the case under consideration the right eigenvectors of Ti coincide with the left 
eigenvectors, \dn) = \di) = \d). However, to keep a general form of the equations, below we will use the right and left 



eigenvectors. In the resonance approximation, the amplitude of the process is defined by equation (3.112), where the 



operators i? 7/ and R 7i are given by equations (|3.116j), E a = e a + (a\T,s E (e a ) + Uyp\a) is the ground state energy 



including the QED corrections of first order in a, Ed and — r^/2 are the real and imaginary parts of an eigenvalue of 

n. 

For the differential cross section in the resonance approximation, we obtain 
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da = (2ir) 4 5(k° f ~ fc?) [ 
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Koleo^A^ldaXdileo^lo)]* 



-2Re 
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d<d' 



(E a + fc° - £ d )2 
o|ea^*|d a )(d£|ea p i^|o)| 2 



r^/4 



£ a + fc? - s d + l r d /2 

a|ea^* |d' fl )* (d' L \ea p A P \a) 



E a + k° i -E d ,-iT d ,/2 
For the total cross section, using the optical theorem, we find 



dkt 



(3.120) 



S 

o-tat = 2(2^) 3 J2 [ 



■Be({a\ea v Air\d R )(d L \ea p A?\a))(T d /2) 
(E a + fc° - i? d )2 + T2/4 



Im((a|ea 1/ AJ"|di i )<dL|ea p Af |o))(£; d - £„ - fc?) 



(£ a + k9 - E d f + T2/4 



(3.121) 



Let us discuss, for simplicity, the c ase s — 2. Only if the d levels have the same quantum numbers, the second term 
on the right-hand side of equation (3.121) is not equal to zero. In the opposite case, which takes place in the process 
under consideration, |d#) = \d^) and, therefore, 



hn.({a\ea v A v i *\d){d\ea p A p i \a)) = 0. 



(3.122) 



It follows that the the total cross section given by eq uation ( 3.121 ) is the sum of Lorentz-type terms. As to the 
differential cross section, the second term in equation (3.120) is not equal to zero even if the states d = 1,2 have 
different quantum numbers. 

Levels close to each other with identical quantum numbers can appear among doubly excited states of high-Z few- 
electron atoms [ p5[ . As an example, we can consider the (2s, 2s)o and (2px/2i ^P\/2)o states of a heliumlike ion which 
can arise in the process of recombination of an electron with a hydrogenlike ion. A detailed theory of this process was 
given in Ref. p7y . The related numerical calculations were presented in Refs. pj£q{. The results obtained in these 
papers are discussed in section IV(D3). 



IV. NUMERICAL EVALUATIONS OF QED AND INTERELECTRONIC-INTERACTION 

CORRECTIONS IN HEAVY IONS 

A. Methods of numerical evaluations and renormalization procedure 



In the previous sections we have demonstrated how formulas for the energy shifts and the transition and scattering 
amplitudes can be derived from the first principles of QED. These formulas usually contain infinite summations over 
intermediate electron states (summations over the bound states and integrations over the continuum) . These sums 
are generally evaluated by using analytical expressions for the Dirac-Coulomb Green function pLE3Lp7 91 or by using 
relativistic finite basis set methods [92 j97f]. In some cases the summation can be performed analytically by employing 
the generalized virial relations for the Dirac equation |)8J . 

Calculations of most QED corrections require the application of a renormalization procedure. To first order in a, 
one has to renormalize the self-energy and vacuum-polarization diagrams (Fig s. 12, 13). 



The renormalized expression for the SE correction is given by equation ( 2.102] ) which implies using the same 
covariant regularization for both terms on the right-hand side. For the numerical evaluation of this correction, it is 
convenient to analytically isolate the ultraviolet divergence in the (a|E(e a )|a) term and cancel it by the counterterm. 
This can be performed by expanding the Dirac-Coulomb Green function in terms of the free Dirac Green function, 



H{\ 



iO)]- 1 = [u- F (l - iO)]- 1 + [w - ff (l - *0)]-V c [w - F (l - iO)]- 1 
+[u> - #0(1 - i0)}- l V c [u> - H(l - ^0)]- 1 Vb[co - H (l - ^0)]- 1 . 



(4.1) 



where Hq = a ■ p + (3m is the free Dirac Hamiltonian. The three terms in equation ([0) inserted into (a|E(e a )|a) 
divide the SE correction into zero-potential, one-potential, and many-potential terms, respectively. The ultraviolet 
dive rgen ces in the zero- and one-potential terms and in the counterterm can be cancelled analytically (see Refs. 
|^9 -101 1 for details). As to the many-potential term, it can easily be shown that it does not contain any ultraviolet 
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divergences. For an overview of other methods of carrying out mass renormalization in numerical calculations, we 
refer to j|. 

The vacu um-po larization correction is determined by equation ( [2.106 ) with the VP potential defined by (2.107). The 
expression ( 2.107 ) is ultraviolet divergent. The simplest way to reno rmalize this expression is to expand the vacuum 
loop in powers of the Coulomb potential by employing equation (4.1). According to the Furry theorem, contributions 
of diagrams with odd numbers of vertices in a vacuum loop (with free Dir ac pr opagators) are equal to zero. Therefore, 
the first non-zero contribution results from second term in the expansion (4.1). Only this contribution, which is called 
the Uehling term, is ultraviolet divergent. This term becomes finite by charge renormalization. The renormalized 
expression for the Uehling potential is given by 
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(4.2) 



where \e\Zp(r) is the density of the nuclear charge distribution (J p(r)dr = 1). The higher order (in Vc) terms are 
finite and their sum is called the Wichmann-Kroll correction [p7[ . However, the regularization is still required in the 
second non-zer o ter m due to a spurious gauge dependent piece of the light-by-light scattering contribution. As was 
shown in J90|]91 ,102], in the calculation of the vacuum polarization charge density based on the partial wave expansion 
of the Dirac-Coulomb Green function, the spurious term does not contribute if the sum over the angular momentum 
quantum number k is restricted to a finite number of terms (|/c| < K). Thus the Wichmann-Kroll contribution can 
be calculated by summing up the partial differences between the full contribution and the Uehling term. 

The renormalization procedures de scribed above can be adopted for calculations of the self-energy and vacuum- 
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polarization screening diagra ms |39 
the bound-electron g factor p,105- 



m 
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104 1 as well as for the QED corrections to the hyperfine splitting and 



B. Energy levels in heavy ions 



1. Hydrogenlike ions 



The relativistic energies of a hydrogenlike ion are determined by the Dirac equation (2.1). For the point-nucleus 
case, the Dirac equation can be solved analytically and the binding energy is given by 
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(4.3) 



where v = n + y/(j + 1/2) 2 — (aZ) 2 — (j + 1/2), n is the principal quantum number, and j is the total angular 
momentum. To obtain the binding energy to higher accuracy, QED and nuclear effects must be taken into account. 
The finite nuclear size correction is calculated by solving the Dirac equation with the potential of an extended 
nucleus and by taking the differe nce betw een the energies for the extended and p oint nucleus models. This can be 
performed numerically (see, e.g., [112, 113|) or, with a good accuracy, analytically [114] . Wi th a relative accuracy of 
~ 0.2% for Z = 1 — 100, this correction is given by the following approximate formulas [114] 
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(4.4) 
(4.5) 



where 7 = y/l — (aZ) 2 , 

f ls {aZ) = 1.380 - 0.162aZ + 1.612(aZ) 2 , 
f 2s (aZ) = 1.508 + 0.2l5aZ + 1.332(aZ) 2 , 
f 2pi/2 (aZ) = 1.615 + 4.319aZ - 9.152(aZ) 2 + 11.87(aZ) 3 , 

and R is an effective radius of the nuclear charge distribution defined by 
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For the Fermi model of the nuclear charge distribution, 
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it is possible to obtain with a very high precision 
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The expectation values for powers of r for a great variety of nuclear charge distribution models are given in Rcf. [ 1 13 1 . 

Next the QED corrections of first order in a should be taken into account. To this order, the QED correction is 
defined by the self-energy and vacuum-polarization diagrams (Figs. 12,13). The energy shift from th e self- energy 
diagram (Fig. 12) combined with the related mass counterterm (Fig. 14) is determi ned by equation (2.102). The 
self-energy correction for heavy ions was fi rst e valuated by Desederio and Johnson [115] who employed a method 
suggested by Brown, Langer, and Schaefer |l!6| . Later, a more efficient method was developed by Mohr [g9| who 
calculated this correction to a very high accuracy in a wide interval of Z. The method of the potential expansion of the 
boun d-electron propagator for the calculation of the SE correcti on to all orders in aZ was developed by Snyderman 
[ 100 1 and numerically realized first by Blundell and Snyderman 1 117]. A very efficient procedure for the self-energy 
calcu lations which is closely related to the methods of Snyderman and Mohr was developed by Yerokhin and Shabaev 
[ 101 1 . An approach in which the ultraviolet divergences are removed by subtractions in coordinate space was worked 
out by Indelicato and Mohr 1 1 IS ] . The method of the partial- wave renormalization for the calculation o f the first-order 
SE correction was developed by Persson, Lindgren, and Salomonson [119| and by Quiney and Grant J120L To date, 
the most accurate ca lcula tions of the SE correction to all orders in aZ w ere p erformed by Mohr |89 121] and by 
Indelicato and Mohr [122] for the point nucleus case and by Mohr and Soff [ 123 1 for the extended nucleus case. This 
correction was comprehensively tabulated by Bei er and co -workers [124| for finite nuclear radii. The highest accuracy 
for low-Z atoms was gained by Jentschura et al. [ 125| , 126| . 

The VP correction (Fig. 13) is the sum of t he U ehling and Wichmann-Kroll contributions. The first contribution 
can easily be calculated using the expression (4.2) for the Uehling potential. Calculations of the Wichmann-Kroll 
contribution to all orders in aZ were performed first by Soff and Mohr [pi| for the extended nucleus case and by 
Manakov et al. [127] f or th e point nucleus case. A comprehensive tabulation of this correction for extended nuclei 
was pre sented in Rcf. [128|. The most accurate calculations for some specific ions were accomplished by Persson et 
al. p9| . 

The QED corrections of second order in a have not yet been calculated completely. Most VP-VP an d SE -VP 
diagrams can be evaluated by the methods developed for the first-order SE and VP corrections (see p,130| and 
references therein). The most difficult task consists in the evaluation of the SE-SE contribution. The simplest 



part of this contribution, th e loo p-after-loop diagram, was calculated by Mitrushenkov and co-workers 1 1 3 1 1 and by 
Mallampalli and Sapirstein |132 for high-Z ions. Thes e numerical calculations were extended to low-Z atoms by 
Mallampalli and Sapirstein |133 [ and by Yerokhin [134|. Recently they were confirmed by analytical calculations 
of Yerokhi n ]135| . As to the residual SE-SE contribution, a specific part of it was e valuated by Mallampalli and 
Sapirstein ||132[] and an estimate of the complete SE-SE c ontribut ion is in progress [136|. For the current status of the 
corresponding calculations for low-Z atoms, we refer to [137,138]. 

The calculations of the corrections discussed above are based on quantum electrodynamics within the external field 
approximation. It means that in these calculations the nucleus is considered only as a source of the external Coulomb 
field Vc ■ The first step beyond this approxi mation consists in evaluating the nu clear r ecoil correction. This correction 
is given by the sum of the lower-order term ( |2. 194 ) and the higher-order term ( p. 195 ). For the point nucleus case, an 
analytical calculation of the lower-order term employing the virial relations for the Dirac equation yields pfl] 
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(4.11) 



The higher-order term was numerically evaluated to all orders in aZ for point nuclei in Refs. 139.14C]. The corre- 
sponding calculations for ex tend ed nuclei were carried out in Refs. [141,142|. In the case of hydrogen, the highest 
accuracy was gained in Ref. [143], 
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Finally, the nuclear polarization correction should be taken into account. This correction results from diagrams 
describing the interaction of the elect ron with the nucleus where the i ntermediate states of the nucleus are excited. 
It was evaluated by Plunien and Soff [144] and by Nefiodov et al. [145|. 

In Table 1 we present the individual contributions to the ground-state Lamb shi ft in 238 u 91+ . The uncertainty 
of the Dirac binding energy results from the uncertainty of the Rydberg constant ]146[ . As can be deduced from 
the table , the present status of the experimental precision on the ground-state Lamb shift in hydrogenlike uranium 
[147-149J provides a test of QED in first order in a on the level of about 5%. 



2. Heliumlike ions 



In heavy heliumlike ions, in addition to the one-electron contributions considered in the previous subsection, the 
two-electron corrections have to be taken into account. To lowest order in a, this correction is defined by the one- 
photon exchange diagram (Fig. 20). The calculation of this diagram causes no problem. To second order in a, we 
should account for the two-photon exchange diagrams (Fig. 21), the self-energy screening diagrams (Fig. 22), and the 
vacuum-polarization screening diagrams (Fig. 23). For the gr ou nd sta te of a hclium likc ion, the contribution of the 
two-photon exchange diagrams is defined by equations ( 2.134| ), ( 2.13S ), and ( [2.140 ). The corresponding expressions 
for the case of a single excited state were obtained in RefipGfl . The derivation of the related formulas for the case of 
degenerate and quasidegenerate states by the TTGF method also c a uses no difficulties. The self-energy and vacuum- 
polarization screening contributions are given by expressions ( 2.153] )-(2.16C) (the renormalization of these expressions 
is considered in detail in Refs. |3i^ , [40| , [42],[43| ) . For the case of quasidegenerate states, the corresponding formulas are 
derived by the TTGF method in Refs. |4jj75(|. 

The two-photon exchange contribution is conveniently divided into two parts. The first part is the one which can 
be derived from the Bre it equation. For the ground state, it is well determined by the lowest-order terms of the 
aZ-expansion series [150], 



A£: (Brcit) = q,2[_ . 15766638 - 0.6356(aZ) 2 ]m . (4.12) 

The second part is the remaining one. Evaluated in Refs. [|7l| , |73| , it was found to be much smaller than the first part. 
The self-energy and vacuum-polarization screening diagrams were evaluated in Refs. p3,BO,103,104|. The related 
calculations for excited states of heliumlike ions were performed for the vacuum-polarization screening diagrams |44| [ 
and, in the case of non-mixed states, for the two-photon exchange diagrams |151|. 

Today, the theoretical uncertainty of the ground-state energy in heavy heliumlike ions is completely defined by the 
uncertainty of the one-electron contribution. In this conne ction , a direct measurement of the two-electron contribution 
to the ground-state energy in heliumlike ions performed in 152 turns out to be rather important. In Table 2 we present 
the individual two-electron contributions to the ground state energy in heliumlike bismuth. The two-photon exchange 
contribution is divided into two parts as described above. The three- and more photon contribution is evaluated 
within the Breit approximation by summing the Z~ l expansion terms for the ground state ene rgy beginning from 
Z~ 3 |39|. For the zerot h ord er in aZ, the coefficients of this expansion are taken from Ref. [153] and for the second 
order in aZ from Ref. [15C ] . The uncertainty of this contribution results from QED corrections yet uncalculated. 
As one can see from the table, to test the second-order QED effects that result from the theory beyond the Breit 
approximation, the experimental precision has to be improved by an order of magnitude. 



3. Lithiumlike ions 



To date, the highest experimental accuracy was reached in Lamb-shift experiments of lithiumlike ions [154-157]. In 
these ions, in addition to the one- and two-electron contributions, the three-electron corrections have to be calculated. 
To second order in a, the three-electron contribution is determined by six diagrams that describe two-photon exchange 
involving all three electrons. One of these diagrams is shown in Fig. 37. In the case of one electron over the closed 
(Is) 2 shell, the unperturbed wave function is a one-determinant function, 
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where a and b denote the core states with opposite signs of the angular momentum projection, v indicates the 
valence state. It can be shown that the derivation of the expressions for the one- and two-electron corrections in 



54 



three-electron atoms is easily reduced to the derivation in one- and two-electron atoms, respectively. As to t he three- 
electron correction of second order in a, it can be derived by the TTGF method using the iden tity ( 2.129 ) and the 
general rules formulated in section II(E5). Such a derivation for the irreducible part yields [158] 
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where P and Q denote the permutations over the outgoing and incoming electron states, respectively; I a bcd(^) = 
(ab\I(u))\cd). The prime at the sum indicates that terms with va nishi ng denominator have to be omitted in the 
summation. The reducible part of the three-electron contribution is |15S] 
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where I a b-cd = I a bcd(sb - £d) - hccd^a. - Sd), A = e v - e a , \x a is the angular momentum projection of the a electron, 
/it, = —fi g, v i s the valence electron with the opposite sign of the angular momentum projection. In derivation of 
equation ( 4.15| ) by the TTGF method, some terms containing the v electron have been cancelled with the corresponding 
terms from the reducible two-electron contribution of the two- photon exchange diagrams (see Ref. [158| for details). 
The accurate QED calculations of all the two- an d thr ee-electron corrections to the 2px/2 — 2s transition energy up to 
second order in a w ere perfo rmed in Refs. p^ . |4^ , [l58|| . Approximate evaluations of these corrections were previously 
considered in Refs. [ 159-164 1 . 

In Table 3 the individual contributions to the 2pi/ 2 — 2s transition energy in lithiumlike uranium are presented. The 
total theoretical value of the transition energy, 280.46(9) ± 0.20 eV, is in agreement with the related experimental 
value, 280.59(10) eV [154|. As can be seen from the table, the first-order QED contribution is -42.93 eV while the 
total second-order QED contribution beyond the Breit approximation amounts to 1.33 ± 0.20 eV. Comparing these 
values with the total theoretical and experimental uncertainties indicates that the present status of the theory for 
lithiumlike uranium provides a test of the QED effects of first order in a on the level of about 0.5% and of the QED 
effects of second order in a, which result from the theory beyond the Breit approximation, on the level of about 15%. 



C. Hyperfine splitting and bound-electron g factor 



1. Hyperfine splitting in hydrogenlike ions 



The ground-state hyperfine splitting of a hydrogenlike ion is conveniently written as [166] 
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Here m p is the proton mass, [i is the nuclear magnetic moment, \xn is the nuclear magneton, and I is the nuclear 
spin. A(aZ) denotes the relativistic factor 
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where 7 = •v/l — (aZ) 2 . S is the nuclear charge distribution correction, e is the nuclear magnetization distribution 
correction (the Bohr-Wcisskopf correction), and £ ra d is the QED correctio n. Th e f ormul as for the first-order QED 
corrections to the hyperfine splitting are derived in the same way as f ormul as ( 2.115 ), (2.117) for the (^-SE corrections. 
For instance, the SE correction is simply determined by equations ( 2.115 ), ( 2.117 ), if 5V is replaced by the hyperfine 
interaction operator 
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and the electron wave functions are replaced by the wave functions of the whole (electron plus nucleus) atomic 
system. The most complete calculations of the QED and nuclear corrections to the hyperfine splitting were presented 
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in [107,110] (see also a recent review in ]167[ ). Table 4 shows the individual contributions to the hyper fine splitt ing in 
hydrogenlike ions. The total theoretical values are compared with the experimental results obtained in 1 168 -171]. The 
uncertainty of the theoretical predictions is mainly determin ed by the uncertainty of the Bohr-Weisskopf correction 
which was evaluated within the single-particle nuclear model [107]. This uncertainty should be considered only as an 
estim ate of the order of magnitude of the real erro r. Except for Ho, the nuclear mag netic moments are taken from 
Ref. |l72| . For Ho the value recommended in 173 is used. In case of 207 Pb, in Ref. |l72| two values of the nuclear 
magnetic mo ment are given. One (/i = 0.592583(9)/ijv) was measured by the nuclear magnetic reso nanc e (NMR) 
method [1741 while another (/i = 0.58219(2)/zjv) results from an optical pumping (OP) experiment [175]. As was 
shown in [176|, the OP value turns out to be very close to that obtained by NMR if it is corrected for an atomic effect 
(see the related discussion in ||). Therefore, in Table 4 the NMR value for the nuclear magnetic moment of lead is 
used. 

Taking into account that the theoretical uncertainties indicated in Table 4 should be considered only as an order 
of magnitude of the real errors, one can deduce that the total theoretical values are in reasonable agreement with the 
experimental ones. However, remeasurements of the nuclear magnetic moments by employing modern experimental 
technique and calculations of the Bohr-Weisskopf effect within many-particle nuclear models are required to promote 
investigations of the hyperfine splitting in hydrogenlike ions. 



2. Hyperfine splitting in lithiumlike ions 

The energy difference between the ground-state hyperfine splitting components of a lithiumlike ion is conveniently 
written as |38|| 
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Here A^- 2s \aZ) denotes the one-electron relativistic factor for the 2s state, 
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S^ 2 ^ is the one-electron nuclear charge distribution correction, £^ 2s ) is the one-electron nuclear magnetization distri- 
bution correction, and x\ j is the one-electron QED correction. The terms B(aZ)/Z and C(aZ)/Z 2 describe the 
interelectronic-interaction corrections of first and second orders in 1/Z, respectively. The first-order interelectronic- 
interaction correction is conveniently derived using the TTGF method with the closed (Is) 2 shell regarded as belonging 
to the vacuum (see section II(E3)). Such a derivation yields f3q] 
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where F and Mp are the total angular momentum of the atom and its projection, v and v' are the valence states 
of electron with quantum numbers (jm) and (j'm r ), respectively; C IMi F , m is the Clebsch-Gordan coefficient, XlMi 
is the nuclear wave function, c and d denote the core states, /i c indicates the angular momentum projection of the 
core electron, and A a {, = e a — £&• Calculations of the nuclear, QED, and i nterelectronic- inte raction corrections to 
the hyperfme splitting in heavy lithiumlike ions were performed in Refs. |3|, 109, 167, 177 -182]. As for hydrogenlike 
ions, the uncertainty of the theoretical values is mainly determined by the un certainty of the Bohr-Weisskopf effect 
evaluated within the single particle nuclear model. However, in Refs. [ 109 , 177 1 it was found that this uncertainty can 
be considerably reduced if the experimental value of the Is hyperfme splitting in the corresponding hydrogenlike ion 
is known. The basic idea of this method is the following. It can be shown that, with a good precision, the ratio of 
the 2s— Bohr-Weisskopf correction to the Is— Bohr-Weisskopf correction is a function of the atomic structure only and 
does not depend on the nuclear structure, 
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For Z=83, f{aZ)-- 
equation 



= 1.078 and, therefore, e^ 2 ^ — 1.078 e^ ls \ The Is— Bohr-Weisskopf correction can be found by the 
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where AE 
A^E D is the l.s 



Dirac 1S ^ ne relativistic value of the Is hyperfme splitting including the nuclear charge distribution correction, 
QED correction, and A_Eexp is the experimental value of the Is hyperfme splitting. For Z=83, this 
method predicts the ground-state hyperfme splitting in lithiumlike bismuth to be 0.7971(2) eV ]179| ] (for comparison, 
the direct evaluation base d on the single-particle nuclear model gives 0.800(4) eV). Recently, this value was confirmed 
by Sapirstein an d Ch eng |183| who obtained 0.79715(13) eV. Both theoretical values agree with the experimental one 
of 0.820(26) eV fl5lj. 



3. Bound-electron g factor 



The bound-electron g factor in a hydrogenlike ion is defined by 
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where fi^ e ' is the operator of the magnetic moment of electron, /i-g is the Bohr magneton, J is the total angular 
momentum of the electron , an d Mj is its projection. For the ground state, a simple relativistic calculation based on 
the Dirac equation yields [184 



.90 



(4/3X1-V1-M) 2 ) 



The QED and nuclear effects give some corrections to this value: 
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Calculation of the nuclear size correction causes no difficulties. For low-Z atoms, this correction can be found by a 
simple analytical f ormula Jl85| | . The QED correction of first order in a was evaluated without expansion in aZ in 
Refs. r|l06|Jl08|JTTlll (see alsoff). 

Direct measurements of the boun d-electro n g factor in hydrogenlike ions are presently being performed by a GSI 
- Univ ersit at Mainz collaboration [186-189|. To date, the experimental result obtained for hydrog enlik e carbon 
(C 5+ ) |189] amounts to g exp = 2.001041596(5) and agrees with the theoretical predictions of Refs. |j,lll|, gtheo = 
2.001041591(7), and of Ref. ]l9C|| , g thco = 2.001041590(2). The collaboration plans to extend these measurements to 
heavy ions. 



Another possibility for investigations of the bound-electron g factor was recently proposed in |191|. In this work it 
was shown that the transition probability between the ground-state hyperfme splitting components of a hydrogenlike 
ion, including the first-order QED and nuclear corrections, is given by 
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where to is the transition frequency, g( e > is the bound-ele ctron g factor defined above, and g™ is the nuclear g factor 
(both g factors are defined to be positive). Formula ( 4.27 ) allows a simple calculation of the QED and nu clear 
corrections to the transition probability using the corresponding corrections to the bound-electron g factor. In [191] 
it was found that in the experimentally int erest ing cases of Pb and Bi, the QED and nuclear corrections increase the 
transition probability by about 0.3% 



In 



192] the lifetime of the upper hypcrfine splitting component in 209 Bi 



209t3;82+ 



was measured to be T cxp =397.5(1.5) /is. Th is result is in good agreement with the theoretical predictions of Rcf. 



[ 191 1 , Ttheo = 399.01 (19)/xs, and of Ref. |193|, T t hco = 398.89/is. Usi ng form ula (4.27) and the experimental values 
of the hyperfine splitting and the transition probability in 209 Bi 82+ [168 192 1 , the experimental value of the bound- 
electron g factor in 209 Bi 82+ is found to be 1.7343(33). The corresponding theoretical value is 1.7310. The individual 
contributions to the bound-electron g factor in 209 Bi 82+ are given in Table 5. ^From this table, it is clear that the 
QED correction has to be included in order to obtain agreement between theory and experiment. 



D. Radiative recombination of an electron with a heavy ion 



In an energetic collision between a high-Z ion and a low-Z target atom, an electron may be captured by the 
projectile, while a simultaneously emitted photon carries away the excess energy and momentum. This process is 
denoted as radiative electron capture (REC). Since a loosely bound target electron can be considered as quasi- free, 
this process is essentially equivalent to radiative recombination (RR) or its time-reversed analogon, the photoelectric 
effect. The relativistic theory of REC was con sidered in detail in [32,194 195|, and the results of this theory are in 
excellent agreement with experiments [196,197]. 

A systematic QED theory of the RR process is described in detail in sections III(D,E) of the present paper. In 
Refs. [ [l4p7]|8(| ] i this theory was employed to study the process of resonance recombination of an electron with a 
heavy hydrogenlike ion in the case of resonance with doubly excited (2s, 2s)o, (2pi/2, 2pi/2)o, (2s, 2p!/2)o,i states of 
the corresponding hcliumlike ion. Later, this theory was used to evaluate QED corrections to radiative recombination 
of an electron with a bare nucleus |}45[ and interelectronic-interaction corrections to radiative recombination of an 
electron with a heavy heliumlike ion |46|. The results of these investigations are briefly discussed below. More details 
can be found in the original papers |34| , |37l , [l5| . [l6| , p6[ . 



1. QED corrections to the radiative recombination of an electron with a bare nucleus 



We consider the radiative recombination of an electron with momentum pi and polarization /i; with a bare nucleus 
that is placed at the origin of the coordinate frame. This corresponds to the projectile system if we study the radiative 
recombination of a free target electron with a bare heavy projectile. To zeroth order in a, the cross section is 
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where 
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(a\ea u A* Lu \p) = (a\ea ■ A* f \p) , 



(4.29) 



\p) = \pi, Hi) is the wave function of the incident electron defined by equation ( 3.6S ), pi — (p°, pi), p\ = y/pf + m 2 is 
the energy of the incident electron, a is the final state of the one-electron atom, kf = (/c?,kj) with k° f and kj being 
the photon energy and momentum, respectively, Vi is the velocity of the incident electron in the nucleus frame. 

The QED corrections of first order in a are defined by diagrams similar to those shown in Fig. 32. The direct 
calculation by the TTGF method yields for the self-energy correction to the amplitude of the process 
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where the mass counterterm has been added and 
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A similar calculation of the VP correction gives 
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Here t/vp(x) is the vacuum-polarization potential defined by equation (2.107) and 

Q"(k a ,z)=-e 2 Jdxdyi, a ( X )Y^ M+) ( X )D pa (k a ,x-y) 
x— / ^Tr[ 7 ^,y,z)7^ + fc ,z,y)]. 
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In addition to these corrections, we have to take into account a contribution originating from changing the photon 
energy in the zeroth-order cross section (4.28) due to the QED correction to the energy of the bound state a. It 
follows that the total QED correction of first order in a to the cross section is given by 
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Here Tqg D = Tg E ; + r V p is the QED correction given by equations fl4.30| ) and ( |4.32j ) . E a and s a are the energies of 
the bound state a with and w ithou t the QED correction, respectively. 

The expressions ( 4.30 ) and ( 4.32 ) contain ultraviolet and infrared divergences. While the ultraviolet divergences can 
be eliminated by the standard renormalization procedure (see the related discussion in section III(B2)), the infrared 
divergences are more difficult to remove. 

The infrared-divergent part of r results from the region of small momenta of the virtual photon and is regularized 
by a non-zero photon mass /i. An evaluation of this part yields 
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The related contribution to the cross section is 



da-„ 



der(°) 



a 



dflf d£lf n 



-21og(/x/m) - log (n/m)* 1 



log 



y/p? + m 2 - jpj 
■y/p? + m 2 + |pi 



(4.35) 



(4.36) 



where ^ is the cross section in the zeroth-order approximation defined by equation (4.28). To cancel the infrared 



divergent contribution (4.36), we have to take into account that any experiment has a finite energy resolution AE. 
It means that any numbers of photons of the total energy less than AE can be emitted in the process. It follows 
that to find the total cross section in the order under consideration, we must include diagrams in which one photon 
of energy k° = ^/k 2 + /i 2 < AE is emitted along with the emission of the photon with the energy k°f ss p\ — e a (we 
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assume AE <C fe2). These diagrams are shown in Fig. 38. Assuming that the energy resolution is sufficiently high 

(AE <C fc2,m), we retain only those contributions from the diagrams shown in Fig. 38 which dominate at AE — > 0. 
Using the standard technique and omitting terms which approch zero at \x — > 0, we obtain for this contribution 
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We can see that the infrared divergent parts in the equations ( 4.36 ) and ( 4.37 ) cancel each other, 



d(7i n f r da 1 (AE) 



dQi 



dQ, 



da^ a 

■K 



dn 



f 



1 - 2 log 2 - 2 log (AE/m) - 



^-F(| Pi |/p°) 



- + log (AE/m) 




(4.39) 



According to this equation, at a fixed energy of the incident electron, the QED correction depends on the photon- 
energy resolution AE and becomes infinite when AE — ► 0. It means that the validity of this equation is restricted 
by the condition — | log (AE/m) | <C 1 . For extension of the t heory beyond this limit it is necessary to include the 



radiative corrections of higher orders in a (see, e.g., 
corrections and removes the singula rity f or AE — > 0. 



It results in an "exponentiation" of the radiative 



In the derivation of the formulas ( 4.37 ), ( 4.39 ) it has been assumed that the incident electrons have a fixed energy. 
These formulas remain also valid in the case when the energy spread of the incident electrons is much smaller than 
the energy i nterval A E in which the photons are detected. However, this is not the case for the present REC 
experiments [196 197 1, where the energy spread of a quasi- free target electron is much larger than the finite photon- 
energy resolution. In that case, the QED correction to the total RR cross section depends on the form of the energy 
distribution of the target electron. Since the form of this distribution is not well determined, the only way to study 
the QED effects in REC processes is to investigate the cross section into a photon-energy interval which is chosen to 
be much larger than the effective energy spread of the quasi-free target electrons and much smaller than the energy 
of the emitted photon. 

The Uehling part of the VP correction to the RR cross section and a part of the SE correction were numerically 
evaluated in Ref. J4J| . Expressed in terms of the unperturbed cross section, the individual QED corrections to the 
total cross section for the radiative recombination into the K-shell of bare uranium are presented in the Table 6. The 
correction <7on results from changing the bound-state energy. It is determined by the term in the square brackets 
of equation (4.34). The correction a h ^ corresponds to the irreducible part of the diagrams describing the first-order 

QED effect on the bound-state electron wave function. The correction <Tc W results from the diagrams describing the 
QED effect on the continuum-state wave function. 

As in bound-state QED, the Uehling approximation may be expected to account for a dominant part of the VP 
correction. As to the self-energy correction, we expect that the terms calculated in Ref. Eq] give a reasonable estimate 
of the order of magnitude of the self-energy correction which is beyond the correction depending on the photon-energy 
interval AE (see equation (4.39)). The relative value of the last correction, which we denote by S(AE), is defined by 
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The photon-energy interval has to be chosen in the range T <C AE <C fc° 
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where T characterizes 



where — Vi/c. 

the energy spread of the incident electrons. 

In the REC experiments which are performed at GSI, the effective electron-energy spread is determined by the 
momentum distribution of the quasi-free target electrons. The width of this spread in the projectile system increases 
with increasing impact energy. In the case of a N2 gas target, which is presently being employed in the experiments, 



(>() 



the effective energy spread in the projectile (heavy ion) frame amounts to about 10-40 keV for the impact energy in 
the range 100-1000 MeV/u. This value will be considerably reduced in the experiments with a H2 gas target which 
are under preparation. In the case of a H2 gas target, to satisfy the conditions on AE given above we can choose 
AE to be 50 keV in the projectile frame for the impact energy 1 GeV/u. The corresponding photon-energy interval 
in the laboratory (gas-target) frame is determined according to the Lorentz transformation, 



AE proi = A£ lab (l - /3 cos lah )/y/l-/3 2 . (4.41) 

At a fixed A£"p ro j, AE\ a b as a function of the polar angle can be found from this equation. For the photon-energy 
interval chosen above, 5(AE) amounts to —0.59% for an impact energy of 1 GeV/u. Adding the Uehling correction 
and the part of the SE correction presented in Table 6 to the correction 6(AE), we find the QED correction to the 
total cross section amounts to —0.92% for an impact energy of 1 GeV/u. For a more accurate evaluation of this effect, 
complete calculations of all the SE corrections are required. 

The results of the numerical evaluation of the differential cross section can be found in Ref. [Q. Here, we note 
only that the differential cross section at the backward direction vanishes at an impact energy close to 130 MeV/u. 
In particular, it results in a relatively large contribution of the QED correction to the backward cross section at an 
energy of 130 MeV/u. At this energy, the QED Gn+bw+AS correction is about 0.022 mbarn/sr while the zeroth-order 
cross section amounts only to 0.009 mbarn/sr. 

2. Interelectronic-mteraction effect on the radiative recombination of an electron with a heavy heliumlike ion 

We consider the non-resonant radiative recombination of an electron with momentum p^ and polarization /Zj with 
a heavy hcliumlike ion in the ground state which is placed at the origin of the coordinate frame. The final state of 
the system is a lithiumlike ion in the state (ls) 2 v, where v denotes the valence state. This picture corresponds to the 
projectile system if we study the radiative recombination of a free target electron with a heavy heliumlikc projectile. 

To zeroth order, the amplitude of the process is given by 

r(°> = (v\ea ■ A)\p) , (4.42) 

where \v) denotes the wave function of the valence electron. The interelectronic- interaction correction of first order in 
1/Z can easily be derived using the TTGF method with the closed (Is) 2 shell regarded as belonging to the vacuum. 
Such a derivation yields M& 
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Here \v) and \c) are the valence and core states, respectively, p® = ypf + m 2 is the energy of the incident electron, 
k°c = pf — e v is the energy of the emitted photon, and fx c indicates the angular momentum projection of the core 



electron. The expressions (4.44) - (4.47) represent the interelectronic- interaction corrections to the amplitude of the 



process. The corresponding corrections to the differential cross section are 
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In addition to this, a contribution originating from a modification of the energy of the emitted photon in the zeroth- 
order cross section due to the interelectronic interaction should be taken into account which is given by 
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Here E v — e v + AS int is the energy of the valence electron including the first-order interelectronic-interaction correc- 
tion, 

A 4!t = EE(- 1 ) p ( PvPc i / ( A ^)i^) . (4.50) 
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The total interelectronic-interaction correction to the cross section in first order in 1/Z is 
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The direct part of the corrections r{ nt and r™* can be accounted for by a modification of the incoming and outcoming 
electron wave functions by the screening potential 

Vc(x)^V c (x)+V 8CT {x), (4.52) 



where 

Kcr(x) = 2a {\[ d vv 2 (aUv) + fliv)) + £ dyy( 9 2 u (y) + ft(y))] ■ (4-53) 

Here g\ s and f\ s are the upper and the lower components of the radial wave function of the ground state, respectively. 
The screening-potential approximation allows one to account for the dominant part of the interelectronic-interaction 
effect and is widely used in practical calculations. 

Numerical results for the interelectronic-interaction correction to the total cross section of radiative recombination 
of an electron with heliumlike uranium are presented in Table 7 (for more extensive data we refer to pq|). The 
calculations are carried out for a capture into the 2s, 2pi/ 2l and 2p 3 / 2 states of lithiumlike uranium and for various 
projectile energies. The results of the rigorous relativistic treatment are compared with the calculations based on 
the screening-potential approximation. A deviation of the complete relativistic results from the screening-potential 
approximation is mainly determined by the term a™ 1 which strongly increases when p® comes close to the resonance 
condition (p° — (e v — e c ) w e n ). Numerical results for the differential cross section are given in Ref. pq| . 

3. Resonance recombination of an electron with a heavy hydrogenlike ion 

We consider the process of recombination of an electron with a very heavy (Z ~ 70 — 110) hydrogenlike ion in 
its ground state for the case of resonance with doubly excited (2s 2s)o, (2j»i/ 2 2p 1 / 2 )o, (2s 2p 1 / 2 )o,i states of the 
corresponding hcliumlikc ion. We assume that the experimentally measured quantity is a part of the total cross 
section which corresponds to the emission of photons with an energy ui rs Ed — E r , where Ed (d = 1, 2, 3, 4) are the 
energies of the doubly excited (2s 2s) , (2pi/ 2 2p 1 / 2 ) , (2s 2p 1 / 2 )o,i states, respectively, and E r (r = 1,2,3,4) are the 
energies of singly excited (Is 2p 1 / 2 )o,i, (Is 2s)o,i states. The amplitude of this process is given by the sum of the 
amplitudes of the dielectronic recombination (DR) and radiative recombination (RR) processes 

Here p° is the energy of the incident electron. Among the doubly excited states {d}, there are states with identical 
quantum numbers ( (2s 2s)o, (2pi/ 2 2p 1 / 2 )o ) , while all the singly excited states {r} have different quantum numbers. 
Main channels of the decay of (Is 2s)i, (Is 2p 1 / 2 )i and (Is 2s)o, (Is 2p 1 / 2 ) states are one- and two-photon transitions 
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to the ground state, respectively. Therefore, two- and three-photon processes give the dominant contribution to the 
cross section. The cross section of the process can be derived using the method described in section III(F). In the 
resonance approximation, such a derivation yields (see Ref. |37j for details) 
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where \ij) = \JMnj'l'p®jl), (n f V) = (1 1/2 0) are the quantum numbers of the Is state of the hydrogenlike ion, 
1= I(e 2s -sis), 
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The cross section (4.54) consists of four terms, a — X)i=i (7 '- The first two terms (<7i, 02) correspond to the DR 



process, the third term (0-3) corresponds to the RR process, and 04 describes the interference between the DR and 
RR processes. The term <72 is caused by the interference of the DR amplitudes on the levels with identical quantum 
numbers (d, d' — 1,2). The magnitude of this term is determined by the overlap of the levels d, d' and can be 
characterized by the nonortogonality integral (d' R \d R ) which is connected with W dd > by the identity 

W dd , = i{£ d -£* d ,){d' R \d R ). (4.60) 



Formula (4.54) was used in 134] for the numerical calculation of the cross section of the resonance recombin ation 



of an electron with hydrogenlike uranium. Later, a more accurate calculation was performed by Yerokhin [ 200 1 . 
According to these calculations, the ratio 02/0 which characterizes the overlap effect amounts up to 30% in the region 
between the maxima of the curve cr(p°). For the parameters of the overlapping levels it was found |34| ] 

\{d R \d' R )\ =0.180, \W dd ,\/\E d - E d ,\ = 0.183. (4.61) 



Using ( 4.60 ) and the identity 

W dl = i(i\(I - P)\d R ) , (4.62) 



the expression (4.54) can be transformed to the following (see Ref. [B7j for details) 
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This expression also consists of four terms a = J2i=i &i- But, in contrast to (4.54), here the Breit-Wigner part of the 
cross section is completely contained in the first term (W±) which is a sum of Lorentz-type terms. The terms <72, CX4 
again correspond to the interference of the DR amplitudes on the levels d, d! — 1,2 and the interference of the DR 
and RR processes, respectively. But, unlike 02, &4, the terms 02, &4 do not contain any admixture of Lorentz-type 
terms. They are given by sums of terms which are odd functions of (Ei s + p® — E d ). (We should note that the 
admixture of Lorentz-type terms in 02, 04 is very small and hardly affects the values 02/0 and 04/0.) The terms 02, 
04 (02, ~&a) lead to a deviation of the shape of the individual resonances from the Lorentz shape. This deviation 
can be characterized by the Low parameter |83J] 

a(E d -E ls + ^)-a{E d -E ls - T -f) 

<j(E d - E ls ) ■ [ -° ' 



Calculation of this parameter for d = 2 in the case of recombination with U 91+ using the results of Refs. [34 200 
yields S « —0.17. The contributions to 6 from the terms W2 and 04 are equal —0.24 and 0.07, respectively. 



Summing the interference terms in (4.63), we obtain a compact formula for 
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The formulas (4.63|) and (4.65|) are more convenient for the numerical calculations than formula (4.54), since they do 



not require any calculation of the radiative amplitudes. These formulas were employed in Rcf. 86jj to calculate the 
resonance recombination of an electron with hydrogenlike lead. 

V. CONCLUSION 

In the present paper we have considered in detail the two-time Green function method for high-Z few-electron 
systems. This method allows one to formulate a perturbation theory for calculations of various physical quantities in 
a rigorous and systematic way in the framework of quantum electrodynamics. To demonstrate the efficiency of the 
method, we have derived formulas for QED and interelectronic-interaction corrections to the energy levels, transition 
and scattering amplitudes in one-, two-, and three-electron atoms. 

For the last decade, the TTGF method was intesively employed in calculations of QED effects in heavy ions. An 
overview of these calculations was also given in th e present paper. Details of the calculations and other applications 



of the method can be found in Refs. []35 4G 158 1. In particular, in Refs. J42|,[43| the vacuum-polarization and self- 



energy screening corrections to the energi es of lithiumlike ions were calculated. The two-photon exchange diagrams 



for lithiumlike ions were evaluated in Ref. [158]. The second-order two-electron diagrams for quasidegenerate states of 
heliumlike ions are studied in |44||75|] . In Rcf. |3^] , the TTGF method was employed to construct an effective-energy 
operator for a high-Z few-electron atom. In Ref. Eq], this method was used to evaluate the QED corrections to 
the radiative recombination of an electron with a bare nucleus. The interelectronic-interaction corrections to the 
radiative electron capture for a heliumlike ion were considered in Ref. [J46f . In Ref. [J80] , the interelectronic-interaction 
corrections to the transition probabilities in heliumlike ions are derived. 

Concluding, the two-time Green function method provides a uniform and very efficient approach for deriving QED 
and interelectronic-interaction corrections to energy levels, transition probabilities, and cross sections of scattering 
processes in high-Z few-electron atoms. Using an effective potential instead of the Coulomb potential of the nucleus 
allows one to extend this approach to many-electron atoms. 
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APPENDIX A: QED IN THE HEISENBERG REPRESENTATION 

In the Heisenberg representation, the basic equations of quantum electrodynamics in the presence of a classical 
time-independent field A^(x) are 

(i —m — e /l c i (x.))i/j(x) = e /I (x)ip(x) — 5mip(x) , 

UA u {x) =jv{x), (Al) 

where j v (x) — (e/2)[ip(x)^^,ip(x)] is the electron-positron current operator. The state vectors in the Heisenberg 
representation are time-independent 

9 t |$) =0. (A2) 

The physical state vectors have to obey a subsidiary condition 

(ft,A"(aO)< + >|*>=0, (A3) 

where (d v A v (x))^ is the positive-frequency part oid v A"(x) . The Heisenberg operators ijj(x), ip(x), and A v (x) obey 
the same equal-time permutation relations as the corresponding free-field operators. However, in contrast to the 
free fields, the permutation relations for arbitrary times remain unknown. Due to the time-translation invariance, 
Heisenberg operators obey the following transformation equation 

exp (iHt)F(0, x) exp (-iHt) = F(t, x) , (A4) 

where H is the Hamiltonian of the system in the Heisenberg representation. For more details, see {2l],[56 81 . 



APPENDIX B: SINGULARITIES OF THE TWO-TIME GREEN FUNCTION IN A FINITE ORDER OF 

PERTURBATION THEORY 

Let us investigate the singularities of the Green function G(E) in a finite order of perturbation theory. To m-th 
order in e, which corresponds to order m/2 in a, it is given by 
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J-oo ml J 



where the label "con" means that disconnected vacuum- vacuum subdiagrams must be omitted. For simplicity, we 
omit here the mass renormalization counterterm. The presence of this term would not change the consideration given 
below. 

Let us consider the contribution of a diagram of m-th order in e. This diagram is defined by a certain order of 



contractions in equati on ( Bl). T he co ntractions between the electron-positron fields and between the photon fields 
give the propagators ( j2.10| ) and ( 2.11 ), respectively. In this Appendix we will use the following representation for 



these propagators 

(0|7> in (x)Vv(y)|0) = S{x° - y°) £ Vn(x)? n (y) exp [-ie n (x° - y )} 

e„>0 

-%° - x°) J2 M*$ n (y) ^P [-ie„(x° - y )} , (B2) 



e„<0 



(0\TK(x)AUy)\0) = -r I 7^ ex P HV^T^|^|]exp[,k.(x- y) ] _ ^ 



V 



(2n) 3 2 v /k r 

Her e, b y definition, 6(i) = (t+ \t\)/(2t) for t ^ and 9(0) = 1/ 2, and a n on-zero photon mass is introduced. Following 



,201 1, we will use the formalism of time-ordered diagrams [ 202,203) to investigate the singularities of Q^ m / 2 \E) 



Let us consider a certain order of the time variables, 

yl>yl_ 1 >---yl>ti>yl_ i >---yi>t>yi_ i >---yl 
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which defines a time-ordered version of the Feynman diagram. The contribution of the Feynman diagram is the sum 
of all time-ordered versions. Each time-ordered version is conveniently represented by a diagram in which th e ve rtices 



are ordered upwards according to increasing time (see, for example, Fig. 39). According to equations (B2) and 



(B3), each electron propagator contains a sum over the whole electron-energy spectrum and each photon propagator 
contains an integral over the photon momentum. Let us place these sums and integrals in front of the expression for 
the time-ordered version of the Feynman diagram under consideration. Then, an electron line is characterized by an 



electron energy e n and, according to (B2), gives a factor (here we are only interested in time-dependent terms) 

exp[-ie n (y?-yl)] = cxp[-i\e n \(y° - y° k )} for e n > 
and 

- exp [ie n (y° - yg)] = - exp [-i\e n \(y° - j/g)] for e n < , 
where in both cases we consider y® > y®. A photon line gives a factor 

exp[-n/k^T^Q/ t -y2)]> 

where we consider again yf > y®. Each time y® in the diagram is marked by a horizontal dashed line. These lines 
may intersect other electron and photon lines (see Fig. 39). For the point of intersection with an electron line, we 
introduce a factor exp (i|e„|y°) exp (— i|e n |y?) = 1, where e n is the energy of the intersected electron. For the point 
of intersection with a photon line we introduce a factor exp (?-\/k 2 + /i 2 yf) exp (— iyk 2 + /i 2 yf) = 1, where k is the 
momentum of the intersected photon. In addition, we represent the factor exp (iE't') exp (—iEt) as 

exp {iE't') exp (-iEt) = exp [i(E' - E)y°J exp [-i(E' - E)y°J ■ ■ ■ 

x exp [i(E' - E)yl] exp [-i(E' - E)y° k ] 
x exp (iE't') exp (-iEt') exp (iEt') exp (-iEy° t _J 
x exp (iEy^Jexp (-iEyl_ a ) ■ • • exp (iEy^ ) exp {-iEt) . (B4) 

As a result of all these representations, the integral over times at fixed intermediate electron states (n) and photon 
momenta (k) is 



/ dyl / *" dy\ m _ x ■■■ 'dt' / dyl_ x ■■■ ' dt dyl 

J -co J -co J —co J —CO J —co J —CO 


«/ — OO 


x exp [i(E' - E)yl\ exp [i(E - E' - ^ \s n \ - ^J A 2 + V 2 )(yl - 

(m) (m) 


-&.J]- 


x exp [i(E -E'-Y, \tn\ - ^ Vk 2 + M 2 )(^ - ?)} 
(0 (0 




x exp [i(E - Y. I £ «l ~ U ^ + <« 2 )(*' - «?,-J] ' ' ' 
(*') (*') 




x exp [i(E Y l £ "l - U V 7 ^ + M 2 )(^ s ~ *)] 




x exp [i(- ^J |e„| - ^ Vk 2 + M 2 )(t - 1&_J] ' ' ' 
(*) (*) 




x exp [f(- Y 1^1 "LVk 2 + M 2 )(2/° - !&)] • 





(B5) 

(2) (2) 

Here X)(m) l £ «l denotes the sum of the electron energies from the electron lines which are sandwiched between the 
horizontal lines corresponding to the times y° and y® _ . X)(m) \/k 2 + /i 2 denotes the sum of the photon energies 
from the photon lines which are sandwiched between the horizontal lines corresponding to the times y® m and y® x - 
Using the identity 



J — oc 



dx exp (— iax) = , (B6) 

a + iO 



G7 



we easily find 



I m = 2nS(E-E')- 



Efm) \ £ n\ ~ Y,( m ) V k2 + A* : 



(m) 



Em M - Em VkMV 



^-ErtoN-E(toV F +^+*o 



s-E w N-E w VkV' 



iO 



-EmM-Em v& 



•/*" 



E 



(2) 



EmV^ 



•/*' 



(B7) 



A similar calculation for t 1 <t yields an expression which is obtained from (B7) by a replacement E — > — -E in all the 
denominators. 



Because each photon line contracts two vertices, at least m/2 denominators in (B7) have to contain the photon- 
energy terms and therefore do not contribute to the singularities under consideration. It follows that C/( m / 2 ) (E) has 
isolated poles of all orders till m/2 + 1 at the unperturbed positions of the bound state energies. The separation of 
these poles from the related cuts is provided by keeping a non-zero photon mass /i. As to the cuts starting from the 
lower energy levels, they can be turned down. 



APPENDIX C: TWO-TIME GREEN FUNCTION IN TERMS OF THE FOURIER TRANSFORM OF THE 

2JV-TIME GREEN FUNCTION 



To prove the equation (2.56) we have to show that 

2tt 1 



n i /-00 

g{E)8{E-E') = -—J dp ,--- dp° N dPi 



dp 



i N! 



N 



P'° N ) 



xG(KV..yM,...,A) 



(Cl) 



where the coordinate variables are omitted, for brevity. According to the definition of G (see equation (2.12)), equation 
(|Cl|) is equivalent to 



■d P % 



Using the identity 



we obtain 



o_ 1 r<x> 

g(E)S(E-E') = ——J dp"--- dp° N dp'l 

x8{E-p\ p° N )S(E'-p' 1 ° p%) 

/oo 
dx° ■ ■ ■ dx%dxf ■ ■ ■ dx% 
-OO 

x exp (i P ?x'° + ■■■ + ip%x% - iplxl ip° N x° N ) 

x (0|I>(!ci) ■ • ■ ^(a^v-)^(arjv) • ■ • #n)|0) 

= W 2N TJn.f dp "--- dp ° Ndp '° ■ ■ ■ dp '° 

CLXi ' ' ' (XX j\j(%X-\ ' ' ' ClX at 



'N 



1>Vl x 2 ' 



xexp[i(E' -pf- ■ 

x exp [-i(E -pi p° N )x{ - ip\x% 

x (0|T^(xi) • • • ^{x' N )^{x N ) ■ ■ • #ci)|0) 



ip N x N ] 
- ip N x N ] 



1 
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du) exp (iojx) = 6(x) 



(C2) 



(C3) 



(»8 



o,_ i roc 

G{E)5{E - E') = (2vr)- 2 — — - / dx° x ■ ■ ■ dx° N dxf ■ ■ ■ dx% 8(x? - x'%) ■ ■ 

xSixl -»§)-.. S(x1 - x° N ) exp (iE'xf - iEx\) 
x (0|T^(xi) • • • *P{x' N Mx N ) ■ ■ • ?(x!)|0) 

= / dx°dx'° exp (iE'x'° - iEx°) 

x (0|TV(z'°, xi) • ■ • if>(xf°, x^(x°, x N ) ■ ■ ■ #r°, xi)|0) 



«-<) 



(C4) 



The last equation exactly coincides with the definition of Q{E) given by ( 2.18|) . 



APPENDIX D: MATRIX ELEMENTS OF THE TWO-TIME GREEN FUNCTION BETWEEN 

ONE-DETERMINANT WAVE FUNCTIONS 



To derive the equation ( p. 61 ) we use the fol lowin g two identities. First, if A is a symmetric operator in the 
coordinates of all electrons, we obtain (see, e.g., [ 204 1) 



A lk = ( Ui \A\u k ) = ^(-l) P Vk(^i) ' ' ' V>k,(£v)^i> • • • »&;&, ■ ■ -,£») 
p 



(Dl) 



where repeated variables {£} imply integration (the integration over x and the summation over a ) and 



Mffii ■ ■ ■ ,CivJ^i) ■ 
in the form 



we can find 



, £jv) is the kernel of the operator A. Second, if the kernel of the operator A is represented 

A(£, . ..,&;&, • • • .60 = E(- 1 ) Qfl ^Qi' • ■ ■ >fe£i, ■ • ■ ,M , (D2) 



A 4fc = N\j2(-i) p r Pl M) ■ ■■r PiN (sM& ■ ■ . ,&;& 
p 

X^fei(Cl)--"0fc N (^Ar)- 



,60 



According to ( 2.6CJ ), we have 

g(B) 7l ° • • • 7 ° n 8(e -e') = ——J dpi--- dpldp? -.-dp 



xS(E-p ( l p%)5{E'-p'° p'») 

x ^(-1) P G((^° 1 , e P1 ), • ■ ■ , {p%, eW); (p?, £i), • • • , (p°n^n)) 



* ■'"• p J-oo 

x«5(£ -p° pawls' - pf p#) 

xG((p' 1 °, & x ), . . . , (p*, £ PJV ); (p?, Ci), . . . , (pSr, 6v)) 
2J(-1) P G(^ P1 , . . . , £ PJV ; £i, . . . , £ w ) . 



Using |5|)-([D|) we easily obtain ( |2.6l| ). 



(D3) 



(D4) 



APPENDIX E: DOUBLE SPECTRAL REPRESENTATION FOR THE TWO-TIME GREEN FUNCTION 

DESCRIBING A TRANSITION PROCESS 



Let us consider the function G(E' , E) defined as 
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/oo 
dtdt' cxp {iE't 1 - iEt)(Q\TA(t')B(0)C(t)\0) . (El) 

-OO 

Using the transformation rules for the Heisenberg operators and integrating over the time variables, we can derive 
the following double spectral representation for G(E' , E) 



/oo 
dW'dW 
-oo 



K(W',W) L(W',W) 



dW'dw 
dwdW 



(E' -W'){E-W) {E 1 + W) (E + W) 

M{W',u>) N(W',oj) 



{E 1 - W')(k° + u>) (E f + W')(k°-u) 
P(w,W) Q(u,W) 



(k° -lj)(E-W) (k° + uj)(E + W) 



(E2) 



where k° = E-E', 

K{W',W) = ^6{W' - E n )6(W - E m )(0\A(0)\n)(n\B(0)\m)(m\C(0)\0) , (E3) 
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L(W, W) =J2 6 ( W ' - E n)HW - E m )(0\C(0)\n)(n\B(0)\m)(m\A(0)\0) , (E4) 

M(W',u) = J2 S( - W ' ~ E ^ 5 ^ ~ E m ){0\A(0)\n){n\C(0)\m){m\B(0)\0) , (E5) 

N(W',u) = J2 S ( W ' - E ™)S(u - £„)<0|B(0)|n)<n|C(0)|m)<m|i4(0)|0) , (E6) 

n,m 

P(«, W) = VJ 6(w - E n )S(W - P m )(0|P(0)|n)(n|^(0)|m)(m|C(0)|0) , (E7) 

Q(w, W) - 51 ^ " ^"O^C^ - P„)(0|C(0)|n)(n|^(0)|m)(m|P(0)|0) . (E8) 
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TABLE I. The ground-state Lamb shift in 238 U 91+ , in eV. 



112 



Point nucleus bindi ng ener gy 
Finite nuclear sizejM, 
First-order SE p3fl ~ 
First-order VP jugf 
Second-order QED 
Nuclear recoil [141 
Nuclear polarization 
Lamb shift theory 



144, 



Lamb shift experiment |14 



-132279.92(1) 

198.81(38) 

355.05 

-88.60 

± 1.5 

0.46 

-0.20(10) 

465.52(39) ± 1.5 

468(13) 



TABLE II. The two-electron contribution to the ground-state energy in 209 Bi 81+ , in eV. 

One-photon exchange contribution 

Two-photon exchange 

within the Breit approximation 

Two-photon exchange 

beyond the Breit approximation 

SE screening 

VP screening 

Three- and more photon contribution 

Total theory |p9| 

Experiment [[152(] 



1897.56(1) 


-10.64 


-0.30(1) 


-6.73 


1.55 


0.06(7) 


1881.50(7) 


1876(14) 



TABLE III. The 2p 1/2 - 2s transition energy in ^ M U by+ , in eV 



One-photon exchange M3| 
One-electron nuclear size 



First-order QED Jll9| , |l23| ] 
Two-photon exchange w ithin 
the Breit appoximation [158 
Two-photon exchange be yond 
the Breit approximation |158| 
Self energy screening p3[ 
Vacuum polarization screening |42l 
Three- and mo re p hoton exchange [165 
Nuclear recoil |139[ 
Nuclear polarization 



144 



One-electron second-order QED 
Total theory 



Experiment [154] 



368.83 

-33.35(6) 

-42.93 

-13.54 

0.17 
1.52 

-0.36 
0.16(7) 
-0.07 
0.03(1) 
±0.20 
280.46(9) ± 0.20 
280.59(10) 



TABLE IV. Theoretical contributions to the ground state hyperfine splitting in hydrogenlike ions (in eV). 



Ion 


l(ibjj 6B+ 


18b Re V4+ 


l« Re 74+ 


aJ7p b si+ 


'2W)-q^+ 


fJ-f/J-N 


4.177(5) 


3.1871(3) 


3.2197(5) 


0.592583(9) 


4.1106(2) 


Relativistic value 












for point nucleus 


2.326(3) 


3.010 


3.041 


1.425 


5.839 


Finite nuclear charge 


-0.106(1) 


-0.213(2) 


-0.215(2) 


-0.149 


-0.649(2) 


Bohr-Weisskopf effect 


-0.020(6) 


-0.034(10) 


-0.035(10) 


-0.053(5) 


-0.061(27) 


QED 


-0.011 


-0.015 


-0.015 


-0.007 


-0.030 


Total theory fjlTSfl 


2.189(7) 


2.748(10) 


2.776(10) 


1.215(5) 


5.100(27) 


Experiment [|l6£-171 


2.1645(5) 


2.719(2) 


2.745(2) 


1.2159(2) 


5.0840(8) 
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TABLE V. The bound-electron g factor in 2UM Bi 



209tj;82+ 



Relativistic value 

QED 

Nuclear size c orre ction 

Total theory jm! 



Experiment [192] 



1.7276 
0.0029 
0.0005 
1.7310 
1.7343(33) 



TABLE VI. The relative values of the QED corrections to the total cross section for the radiative recombination into the 
K-shell of bare uranium [45], expressed in %. 



Impact energy [MeV/u] 


Correction 


Vacuum 


polarization, 


in% 


Self 


energy, in % 


100 


a W 

u CW 

Total 




0.126 

-0.006 

0.120 






-0.390 

? 
-0.390 


300 


a W 

u CW 

Total 




0.175 

-0.003 

0.173 






-0.513 

? 
-0.513 


1000 


a W 

u CW 

Total 




0.220 
0.043 
0.263 






-0.591 

7 

-0.591 



TABLE VII. The zeroth-order cross section a^°' and the first-order interelectronic-interaction correction calculated in Ref. 



[46], in barns. aU denotes the interelectronic-interaction correction calculated in the screening potential approximation and 
°int indicates the results of the rigorous relativistic calculation. 



Impact energy [Mev/u] 



r(°) 



a (1) 

u scr 



r w 



2s-state: 



100 
300 
700 



41.203 
9.105 

2.457 



-1.393 

-0.3345 

-0.0979 



-2.055 

-0.3755 

-0.1051 



2pi/2-state: 



100 
300 
700 



33.041 
5.042 
1.065 



-2.535 

-0.4538 

-0.1022 



-3.088 

-0.3864 

-0.0861 



2p 3 / 2 -state: 



100 
300 
700 



31.489 
3.646 
0.622 



-2.275 

-0.3132 

-0.0568 



-2.896 

-0.2804 

-0.0489 
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FIG. 1. Singularities of the two-time Green function in the complex E plane. 
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FIG. 2. Singularities of the two-time Green function in the bound-state region, if the interaction between the electron-positron 
field and the electromagnetic field is switched off. 
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FIG. 3. Singularities of the two-time Green function in the bound-state region, disregarding the instability of excited states. 
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FIG. 4. Singularities of the two-time Green function in the bound-state region if the cuts are turned down, to the second 
sheet of the Riemann surface. The instability of excited states is disregarded. 
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FIG. 5. Singularities of the two-time Green function in the bound-state region if the cuts are turned down, to the second 
sheet of the Riemann surface. The instability of excited states is taken into account. 
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FIG. 6. Singularities of the two-time Green function in the bound state region for a non-zero photon mass, including one- 
and two-photon spectra, if the cuts are turned down, to the second sheet of the Riemann surface. The instability of excited 
states is disregarded. 
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FIG. 7. The contour V surrounds the pole corresponding to the level under consideration and keeps outside all other singu- 
larities. For simplicity, only one- and two-photon spectra are displayed. 
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FIG. 8. Singularities of the photon propagator in the complex u> plane for a non-zero photon mass /i. 



FIG. 9. The contour V surrounds the poles corresponding to the quasidegenerate levels under consideration and keeps outside 
all other singularities. For simplicity, only one-photon spectra are displayed. 
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FIG. 10. A deformation of the contour F that allows drawing the cuts to the related poles in the case of quasidegenerate 
states when /i — > 0. For simplicity, only one-photon spectra are displayed. 



FIG. 11. The interaction with an external potential 8V(x) 




FIG. 12. The first-order self-energy diagram. 
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FIG. 13. The first-order vacuum-polarization diagram. 
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FIG. 14. The mass-counterterm diagram. 
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FIG. 15. SV - self-energy diagrams. 
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FIG. 16. SV - mass-counterterm diagrams. 
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FIG. 17. C is the contour of the integration over the electron energy u) in the formalism with the standard vacuum. C" is 
the integration contour for the vacuum with the (Is) 2 shell included. The integral along the contour Cint — C — C describes 
the interaction of the valence electron with the (Is) 2 — shell electrons. 





FIG. 18. The first-order self-energy diagrams for a two-electron atom. 
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FIG. 19. The first-order vacuum-polarization diagrams for a two-electron atom. 



J\nj\TL 



FIG. 20. One-photon exchange diagram. 
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FIG. 21. Two-photon exchange diagrams. 
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FIG. 22. The self-energy screening diagrams. 
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FIG. 23. The vacuum-polarization screening diagrams. 



85 



JVlA/l J~\J~U~U~\- 



j-uxnsx. 



JVU1A. •- 



FIG. 24. SV - interelectronic-interaction diagrams. 



FIG. 25. Coulomb nuclear recoil diagram. 



FIG. 26. One-transverse-photon nuclear recoil diagrams. 
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FIG. 27. Two-transverse-photon nuclear recoil diagram. 



FIG. 28. Two-electron Coulomb nuclear recoil diagram. 



FIG. 29. Two-electron one-transverse-photon nuclear recoil diagrams. 
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FIG. 30. Two-electron two-transverse-photon nuclear recoil diagram. 
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FIG. 31. The photon emission by a one-electron atom in zeroth-order approximation. 
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FIG. 32. The first-order QED corrections to the photon emission by a one-electron atom. 
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FIG. 33. The mass-counterterm corrections to the photon emission by a one-electron atom. 
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FIG. 34. The radiative recombination of an electron with a hydrogenlike atom in zeroth-order approximation. 
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FIG. 35. The interelectronic-interaction corrections of first order in 1/Z to the radiative recombination of an electron with 
a hydrogenlike atom. 
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FIG. 36. The photon scattering on a one-electron atom in zeroth-order approximation. 
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FIG. 37. A typical diagram describing two-photon exchange between three electrons in a lithiumlike atom. 
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FIG. 38. The radiative recombination accompanied by emission of a soft photon. 
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FIG. 39. A time-ordered version of a Feynman diagram. 
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